
/*

	Replication code for McGuirk, Hilger, and Miller, "No Kin in the Game: Moral Hazard and War in the US Congress"
	
	Journal of Political Economy, 2023
	
	Email: eoinfmcguirk@gmail.com or eoin.mcguirk@tufts.edu, nathaniel.hilger@gmail.com, nicholas.l.miller@dartmouth.edu
	
	Instructions
		1. Install additional packages if necessary (line 36)
		2. Insert filepath (line 29)
		3. Run file
	
*/


***************************************************
***************************************************
***************************************************
*	PROLOGUE
***************************************************
***************************************************
***************************************************


*	Declare folder [insert filepath]
clear
global folder = "[INSERT FILEPATH HERE]/mhm_nokin_JPE_REP_FILES"
global folder = "/Volumes/GoogleDrive/My Drive/Moral Hazard and War/7_JPE_REPLICATION/mhm_nokin_JPE_REP_FILES"

*	Set directory 
cd "$folder"

*	Set Stata version
version 17

*	Install additional packages if necessary
/*
ssc install reghdfe, 	replace
ssc install rdrobust, 	replace
ssc install rdplot, 	replace
ssc install eventdd, 	replace
ssc install coefplot, 	replace
*/

*	Set output directory
global output =  "$folder/output"

*	Set up 
set more off
set seed 12345
set scheme lean1
log using "$folder/output/log.txt", text

/*	Contents

*	A.	REPLICATION CODE: MAIN VOTES
*	B.	REPLICATION CODE: HAWKS AND DOVES

*/

***************************************************
***************************************************
***************************************************
*	A.	REPLICATION CODE: MAIN VOTES
***************************************************
***************************************************
***************************************************

*	Load data
use "mhm_nokin_JPE_REP_DATA.dta", clear

*	Restrict to main vote data 
keep if vote_pro_cons !=. // 26,373

*--------------
*		TABLE 1
*--------------
		est clear
		local X terms age age2 i.party_vv house0 
		
			*	1.		
			reghdfe vote_pro_cons  	sons_any_draft_running_b4 ///
				if kids_any_draft_running_b4 == 1, a(kids_num_now sons_num_now) cluster(id_vv voteFE)
						
			estadd local lfe "No"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			est sto reg01
 
 			*	2.		
			reghdfe vote_pro_cons  	sons_any_draft_running_b4 	`X'	///
				if kids_any_draft_running_b4 == 1, a(kids_num_now sons_num_now) cluster(id_vv voteFE)
 		
			estadd local lfe "No"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			est sto reg02
			
			*	3.		
			reghdfe vote_pro_cons  	sons_any_draft_running_b4 ///
				if kids_any_draft_running_b4 == 1, a(kids_num_now sons_num_now voteFE) cluster(id_vv voteFE)
						
			estadd local lfe "No"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			est sto reg03
 
 			*	4.		
			reghdfe vote_pro_cons  	sons_any_draft_running_b4 	`X'	///
				if kids_any_draft_running_b4 == 1, a(kids_num_now sons_num_now voteFE) cluster(id_vv voteFE)
 		
			estadd local lfe "No"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			est sto reg04
			
		esttab  reg01 reg02 reg03 reg04 using "$output/JPE_REP_T1.tex", replace se noconstant keep(sons_any_draft_running_b*  ) order(sons_any_draft_running_b*  ) ///
		star(* 0.100 ** 0.050 *** 0.010) b(%10.4f %10.4f)  nonotes   ///
		scalars("con 	Controls"				///
				"vfe 	Vote FE"   				///
				"sfe 	Number of Sons FE"		///
				"kfe 	Number of Children FE"	/// 
				"clust1 Legislators" 			///
				"clust2 Votes" 					///
				"ymean Mean Dep. Var." 			///
				"N Observations")  				///
		mgroups("Pro Draft Vote" , pattern(1 0 0 0)	///
		prefix(\multicolumn{@span}{c}{) suffix(})  	///
		span erepeat(\cmidrule(lr){@span}))  		///
		sfmt(%~12s %~12s %~12s %~12s %~12s %~12s %10.2f %10.0f)  /// 
		label obslast nomtitles booktabs
						
			
*--------------
*		FIG 1
*--------------
		rdrobust vote_pro_con 	son_RD  	if  sons_num_now ==1 ,   vce(cluster id_vv)  masspoints(adjust)
		rdplot 	 vote_pro_con 	son_RD 	if  sons_num_now ==1 & 	inrange(son_RD, - 6,  6)   , h( 6) bw(mserd)  p(1)  c(0)  vce(cluster id_vv) ///
		graph_options(title("") ylabel(0 0.2 0.4 0.6 0.8 1 ) xlabel(-6(1)6)	 	///
		ytitle(Pro Draft Vote) xtitle(Son's Age {&minus} Upper Cutoff) legend(off)) masspoints(adjust)  ci(95)  shade 
		graph save "$output/JPE_REP_F1a", replace 
		graph export "$output/JPE_REP_F1a.png", as(png)   replace
		graph export "$output/JPE_REP_F1a.pdf", as(pdf)   replace
								
		rdrobust vote_pro_con 	dtr_RD  if  dtrs_num_now ==1,   vce(cluster id_vv)  masspoints(adjust)
		rdplot 	 vote_pro_con 	dtr_RD 	if  dtrs_num_now ==1 & 	inrange(dtr_RD, - 6,  6)   , h( 6) bw(mserd)  p(1)  c(0)  vce(cluster id_vv) ///
		graph_options(title("") ylabel(0 0.2 0.4 0.6 0.8 1 ) xlabel(-6(1)6)	 	///
		ytitle(Pro Draft Vote) xtitle(Daughters's Age {&minus} Upper Cutoff) legend(off)) masspoints(adjust)  ci(95) shade
		graph save "$output/JPE_REP_F1b", replace 
		graph export "$output/JPE_REP_F1b.png", as(png)   replace
		graph export "$output/JPE_REP_F1b.pdf", as(pdf)   replace		
					
				
*--------------
*		FIG 2
*--------------			
		est clear
		cap est sto drop

		forvalues i = 4(1)10 {
			rdrobust vote_pro_ 	son_RD  if  sons_num_now ==1, masspoints(adjust)  vce(cluster id_vv) h(`i')
				est sto B`i'
		}

		coefplot  B4  ||B5 ||B6 ||B7 ||B8 ||B9|| B10 ,  yline(0,  lcolor(red))  bycoefs byopts(yrescale) vertical xtitle(Bandwidth) ytitle(RD Estimate) ///
		ylabel(-0.3(0.1)0.5)name("os", replace)  msymbol(circle) ysize(16) xsize(24) ytitle("RD Estimates of {&rho}", margin(0 -8 0 0))			
					
		graph save 		"$output/JPE_REP_F2a", replace 
		graph export	"$output/JPE_REP_F2a.png", as(png)   replace
		graph export 	"$output/JPE_REP_F2a.pdf", as(pdf)   replace			
			
			
		forvalues i = 4(1)10 {
			rdrobust vote_pro_ 	dtr_RD  if  dtrs_num_now ==1, masspoints(adjust)  vce(cluster id_vv) h(`i')
				est sto D`i'
		}

		coefplot   D4 || D5 || D6 || D7 || D8 || D9||  D10 ,  yline(0,  lcolor(red))  bycoefs byopts(yrescale) vertical xtitle(Bandwidth) ytitle(RD Estimate) ///
		ylabel(-0.3(0.1)0.5) name("os", replace)  msymbol(circle) ytitle("RD Estimates (Placebo)", margin(0 -8 0 0))	
					
		graph save 		"$output/JPE_REP_F2b", replace 
		graph export	"$output/JPE_REP_F2b.png", as(png)   replace
		graph export 	"$output/JPE_REP_F2b.pdf", as(pdf)   replace				
			
	
*--------------
*		TABLE 2
*--------------
		est clear
		local X terms age age2 i.party_vv house0 
		
			*	1.
			reghdfe vote_pro_cons  	sons_any_draft_running_b4  kids_any_draft_running_b4 ///
				if (everSTb4  ==1 | everSCb4  ==1), a(kids_num_now sons_num_now id_vv) cluster(id_vv voteFE)
						
			estadd local lfe "Yes"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg01
 
			*	2.
			reghdfe vote_pro_cons  	sons_any_draft_running_b4  kids_any_draft_running_b4 `X' ///
				if (everSTb4  ==1 | everSCb4  ==1), a(kids_num_now sons_num_now id_vv) cluster(id_vv voteFE)
						
			estadd local lfe "Yes"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg02
			
			*	3.
			reghdfe vote_pro_cons  	sons_any_draft_running_b4  kids_any_draft_running_b4 ///
				if (everSTb4  ==1 | everSCb4  ==1), a(kids_num_now sons_num_now id_vv voteFE) cluster(id_vv voteFE)
						
			estadd local lfe "Yes"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg03
 
			*	4.
			reghdfe vote_pro_cons  	sons_any_draft_running_b4  kids_any_draft_running_b4 `X' ///
				if (everSTb4  ==1 | everSCb4  ==1), a(kids_num_now sons_num_now id_vv voteFE) cluster(id_vv voteFE)
						
			estadd local lfe "Yes"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg04
			
		esttab  reg01 reg02 reg03 reg04 using "$output/JPE_REP_T2.tex", replace se noconstant keep(sons_any_draft_running_b* kids_any_draft_running_b*) order(sons_any_draft_running_b* kids_any_draft_running_b*) ///
		star(* 0.100 ** 0.050 *** 0.010) b(%10.4f %10.4f)  nonotes   ///
		scalars("con 	Controls"				///
				"vfe 	Vote FE"   				///
				"lfe 	Legislator FE"			///
				"sfe 	Number of Sons FE"		///
				"kfe 	Number of Children FE"	/// 
				"clust1 Legislators" 			///
				"clust2 Votes" 					///
				"ymean Mean Dep. Var." 			///
				"N Observations")  				///
		mgroups("Pro Draft Vote" , pattern(1 0 0 0)	///
		prefix(\multicolumn{@span}{c}{) suffix(})  	///
		span erepeat(\cmidrule(lr){@span}))  		///
		sfmt(%~12s %~12s %~12s %~12s %~12s %~12s  %~12s  %10.2f %10.0f)  /// 
		label obslast nomtitles booktabs

	 
*--------------
*		FIG 3
*--------------					
		eventdd vote_pro_cons     i.youngest27_kid_event*  , 												timevar(youngest_son_event ) leads(5) lags(5) accum ///
		method(hdfe, absorb(id_vv voteFE  kids_num_now  sons_num_now) cluster(id_vv voteFE)  ) baseline(-1)   level(95) noline   graph_op(ytitle("Event Study Estimates of {&Phi}{superscript:s}{subscript:j}", margin(0 -8 0 0)) xtitle("Son's Age Relative to Upper Cutoff")  ylabel(-0.15(0.05)0.3) xlabel(-5(1)5) legend(off)) endpoints_op(msymbol(circle)) 
		graph save "$output/JPE_REP_F3a", replace 
		graph export "$output/JPE_REP_F3a.png", as(png)   replace
		graph export "$output/JPE_REP_F3a.pdf", as(pdf)   replace

					
		eventdd vote_pro_cons   	age terms age2 i.party_vv i.house0  i.youngest27_kid_event* , 	timevar(youngest_son_event ) leads(5) lags(5) accum ///
		method(hdfe, absorb(id_vv voteFE  kids_num_now  sons_num_now) cluster(id_vv voteFE)  ) baseline(-1)   level(95) noline   endpoints_op(msymbol(circle)) graph_op(ytitle("Event Study Estimates of {&Phi}{superscript:s}{subscript:j}", margin(0 -8 0 0)) xtitle("Son's Age Relative to Upper Cutoff") ylabel(-0.15(0.05)0.3)  xlabel(-5(1)5) legend(off))
		graph save "$output/JPE_REP_F3b", replace 
		graph export "$output/JPE_REP_F3b.png", as(png)   replace
		graph export "$output/JPE_REP_F3b.pdf", as(pdf)   replace
	 	

*--------------
*		FIG 4a -- HAWKS AND DOVES
*--------------	

*--------------
*		FIG 4b & TABLE A17
*--------------	

		est clear
		local X terms age age2 i.party_vv house0 
		
		foreach i in 1 2 5 {
								
			reghdfe win  sons_any_draft_running_b4 `X'	///
				if warFE`i' == 1 & kids_any_draft_running_b4 == 1, a(kids_num_now  sons_num_now) cluster(id_vv voteFE)	
	 		
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			est sto reg`i'1
					

			reghdfe win  sons_any_draft_running_b4 `X'	///
				if warFE`i' == 1 & kids_any_draft_running_b4 == 1, a(kids_num_now  sons_num_now voteFE) cluster(id_vv voteFE)	
	 		
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg`i'2							
	}
	
		esttab  reg11 reg12  reg21 reg22  reg51 reg52  using "$output/JPE_REP_TA17.tex", replace se noconstant keep(sons_any_draft_running_b*) order(sons_any_draft_running_b* ) ///
		star(* 0.100 ** 0.050 *** 0.010) b(%10.4f %10.4f)  nonotes   ///
		scalars(	"vfe \\ Vote FE"  			///
					"sfe Number of Sons FE" 	///
					"kfe Number of Children FE" ///
					"con Other Controls" 		/// 
					"clust1 Legislators"  		///
					"clust2 Votes" 				///
					"ymean Mean Dep. Var." 		///
					"N Observations")  			///
		mgroups("WW1" "WW2" "Cold War" , pattern(1 0  1 0  1  0)  	///
		prefix(\multicolumn{@span}{c}{) suffix(})  					///
		span erepeat(\cmidrule(lr){@span}))  						///
		sfmt(%~12s  %~12s %~12s %~12s  %~12s  %~12s %10.3f %10.0f)  	/// 
		label obslast nomtitles booktabs
	 
		coefplot 	reg11 reg12, keep(sons_any_draft_running_b4)  bylabel("World War I")  || ///
					reg21 reg22, keep(sons_any_draft_running_b4)  bylabel("World War II") || ///
					reg51 reg52, keep(sons_any_draft_running_b4)  bylabel("Cold War") ||,  ///
					byopts(compact cols(1)) xline(0) subtitle(, size(medium) margin(small) justification(left) color(black) bcolor(white) bmargin(top_bottom)) ///
					legend(rows(1)) xtitle("Impact of Draft Age Son on Pr(Reelection)")  levels(90)  xline(0) xlabel(-0.5(0.1)0.2) msymbol(o)

		/*
		graph save  "$output/JPE_REP_F4b.gph", replace
		graph export "$output/JPE_REP_F4b.png", as(png) replace
		graph export "$output/JPE_REP_F4b.pdf", as(pdf) replace
		*/

		
*--------------
*		TABLE A1
*--------------	
		preserve
		*	Main Sample Votes
			use mhm_nokin_JPE_REP_DATA.dta, clear 
			keep if vote_pro !=. 

			gen cf_con_num_draftvotes = . 
			gen cf_pro_draft = . 
			forvalues i = 0(1)1 {
			foreach j in  65 66 76 77 79 82 88 89 90 91 92 93 {
				
			quietly tab votenum if house == `i' & congress == `j'
			replace cf_con_num_draftvotes = r(r) if house == `i' & congress == `j'
			
			quietly sum vote_pro if house == `i' & congress == `j'
			replace cf_pro_draft = r(mean) if house == `i' & congress == `j'
			}
			}
			
		*		Margin
			replace aye = 0 if aye !=1 
			bys congress house0 votenum: egen mean_aye = mean(aye) 
			bys congress house0 votenum: gen gap = abs(-1 + 2*mean_aye)
			bys congress house0: egen cf_con_meangap_draftvotes = mean(gap) 

		*		Save
			bys congress house0: keep if _n == 1
			keep congress house0 cf_con_num_draftvotes cf_con_meangap_draftvotes cf_pro_draft
			save TA1_con_draft_confident_sumstats, replace
		restore


		preserve
		*	All Draft Votes
			use mhm_nokin_JPE_REP_DATA.dta, clear 
			keep if all_draft ==1 

			gen con_num_draftvotes = . 
			forvalues i = 0(1)1 {
			foreach j in  65 66 76 77 79 82 88 89 90 91 92 93 {
			quietly tab votenum if house == `i' & congress == `j'
			replace con_num_draftvotes = r(r) if house == `i' & congress == `j'
			}
			}

		*		Margin
			replace aye = 0 if aye !=1 
			bys congress house0 votenum: egen mean_aye = mean(aye) 
			bys congress house0 votenum: gen gap = abs(-1 + 2*mean_aye)
			bys congress house0: egen con_meangap_draftvotes = mean(gap) 

		*		Save
			bys congress house0: keep if _n == 1
			keep congress house0 con_num_draftvotes con_meangap_draftvotes con_num_draftvotes
			save TA1_con_draft_sumstats, replace
		restore


		preserve
		*	Assemble
			use TA1_con_draft_sumstats, clear

			merge 1:1 congress house0 using TA1_con_draft_confident_sumstats
			drop _merge 

			order cf_con_num_draftvotes con_num_draftvotes   cf_pro_draft cf_con_meangap_draftvotes con_meangap_draftvotes 
			egen chouse = group(congress house)

			label define house0 0 "House", modify
			label define house0 1 "Senate", modify

			drop if congress <= 64
			drop if congress == 78
			replace congress = - congress

			format  cf_pro_draft  con_meangap_draftvotes cf_con_meangap_draftvotes  %10.2f

			latabstat   cf_con_num_draftvotes cf_pro_draft cf_con_meangap_draftvotes   con_num_draftvotes  con_meangap_draftvotes   ///
			if house == 1  , by(congress) tf("$output/JPE_REP_TA1") replace cap("Summary Statistics: Votes") format 
			
			latabstat   cf_con_num_draftvotes cf_pro_draft cf_con_meangap_draftvotes   con_num_draftvotes  con_meangap_draftvotes   ///
			if house == 0  , by(congress) tf("$output/JPE_REP_TA1") append cap("Summary Statistics: Votes") format 
		*		Add totals
					latabstat   cf_con_num_draftvotes cf_pro_draft cf_con_meangap_draftvotes   con_num_draftvotes  con_meangap_draftvotes ///
					if house == 1, s(sum sd) cap("Summary Statistics: Votes") format 
					latabstat   cf_con_num_draftvotes cf_pro_draft cf_con_meangap_draftvotes   con_num_draftvotes   con_meangap_draftvotes ///
					if house == 0, s(sum sd) cap("Summary Statistics: Votes") format 
		restore

		erase TA1_con_draft_confident_sumstats.dta
		erase TA1_con_draft_sumstats.dta
		
*--------------
*		FIG A1
*--------------	
		preserve
		
			bys congress house0 votenum : keep if _n == 1

			label var running_upper "Upper Cutoff"
			label var running_lower "Lower Cutoff"

			sort house congress votenum
			gen vote_time = _n
			tsset vote_time


			tsline  running_upper running_lower if house0 == 0, tit("House of Representatives")xlab(none) ytit("Age") xtit("Votes in Chronological Order") name(JPE_REP_FA1a, replace) ///
				xmla(3 "World War 1"  11 "World War 2" 18 "Korea" 26 "Vietnam", tlcolor(grey) labsize(small) labcolor(grey)) ///
				xli(6, lcolor(grey) lw(vthin)) ///
				xli(17, lcolor(grey) lw(vthin)) ///
				xli(19, lcolor(grey) lw(vthin))	legend(ring(0) position(2) )
				
			graph save "$output/JPE_REP_FA1a", replace 
			*graph export "$output/JPE_REP_FA1a.png", as(png) replace 
			*graph export "$output/JPE_REP_FA1a.pdf", as(pdf) replace 
			
			tsline  running_upper running_lower if house0 == 1, tit("Senate")xlab(none) ytit("Age") xtit("Votes in Chronological Order") name(JPE_REP_FA1b, replace)  ///
				xmla(45 "World War 1"  71 "World War 2" 91 "Korea" 118 "Vietnam", tlcolor(grey) labsize(small) labcolor(grey)) ///
				xli(56, lcolor(grey) lw(vthin)) ///
				xli(87, lcolor(grey) lw(vthin)) ///
				xli(95, lcolor(grey) lw(vthin)) legend(ring(0) position(2) )

			graph save "$output/JPE_REP_FA1b", replace 
			*graph export "$output/JPE_REP_FA1b.png", as(png) replace 
			*graph export "$output/JPE_REP_FA1b.pdf", as(pdf) replace 
			
			graph combine "$output/JPE_REP_FA1a" "$output/JPE_REP_FA1b", cols(1)
			graph save "$output/JPE_REP_FA1", replace 
			graph export "$output/JPE_REP_FA1.png", as(png)   replace
			graph export "$output/JPE_REP_FA1.pdf", as(pdf)   replace
			
		restore

		
*--------------
*		TABLE A2
*--------------	
		estpost tab running_lower, miss
		esttab, cell("b pct(fmt(a))")  collab("Freq." "Percent")  noobs nonumb nomtitle tex

		estpost tab running_upper, miss
		esttab, cell("b pct(fmt(a))")  collab("Freq." "Percent")  noobs nonumb nomtitle tex

						
*--------------
*		TABLE A3
*--------------	

		preserve

			label var sons_any_draft_running_b0 	"\hspace{5mm} Draft Age Son"
			label var dtrs_any_draft_running_b0 	"\hspace{5mm} Draft Age Daughter"
			label var kids_any_draft_running_b0 	"\hspace{5mm} Draft Age Child"
			label var sons_any_draft_running_b4 	"\hspace{5mm} Draft Age Son (4 Year Lead)"
			label var dtrs_any_draft_running_b4 	"\hspace{5mm} Draft Age Daughter (4 Year Lead)"
			label var kids_any_draft_running_b4 	"\hspace{5mm} Draft Age Child  (4 Year Lead)"
			label var sons_any_draft_window_b0 		"\hspace{5mm} Draft Age Son"
			label var dtrs_any_draft_window_b0 		"\hspace{5mm} Draft Age Daughter"
			label var kids_any_draft_window_b0 		"\hspace{5mm} Draft Age Child"
			label var sons_any_draft_window_b2 		"\hspace{5mm} Draft Age Son (2 Year Lead)"
			label var dtrs_any_draft_window_b2 		"\hspace{5mm} Draft Age Daughter (2 Year Lead)"
			label var kids_any_draft_window_b2 		"\hspace{5mm} Draft Age Child (2 Year Lead)"

			* Legislator-Vote Level
			estpost sum vote_pro running_upper_cutoff running_lower_cutoff age terms ///
			dem1 house0 running_vote sons_any_draft_running_b0 dtrs_any_draft_running_b0 ///
			kids_any_draft_running_b0 sons_any_draft_running_b4 dtrs_any_draft_running_b4 ///
			kids_any_draft_running_b4  window_vote sons_any_draft_window_b0 dtrs_any_draft_window_b0 ///
			kids_any_draft_window_b0 sons_any_draft_window_b2 dtrs_any_draft_window_b2 ///
			kids_any_draft_window_b2 if vote_pro_cons !=.  , d
			
			est store ds_LV_main

			esttab ds_LV_main using "$output/JPE_REP_TA3b", var(5) replace ///
						mgroups("Legislator-Vote Level Variables, 140 Votes", pattern(0 1 )  ///
						prefix(\multicolumn{@span}{c}{) suffix(})  ///
						span erepeat(\cmidrule(lr){@span}))  ///
			cells("mean(fmt(2) label(Mean)) sd(fmt(2) label(SD)) count(fmt(%9.0fc) label(Count)) min(fmt(%9.0fc) label(Min)) p50(fmt(%9.0fc) label(Median))  max(fmt(%9.0fc) label(Max))" ) ///
			label  nonumber noobs booktabs
			
		restore 
		
		preserve
		
			label var income 		"Income in 1939 (From 1940 Census)"		
			bys  id_vv: gen tag_L = _n==1
			tab tag_L
			
			estpost sum 	kids_nonmissing_data row_total_missing_yob row_total_missing_gender  		///
							kids_any_total_nonmiss kids_num_total sons_any_total_nonmiss sons_num_total ///
							dtrs_any_total_nonmiss  dtrs_num_total   income  if tag_L == 1 & vote_pro_cons !=.  , d

			est store ds_L_main

			esttab ds_L_main using "$output/JPE_REP_TA3a", var(5) replace 	///
				mgroups("Legislator Level Variables", pattern(0 1 )  		///
				prefix(\multicolumn{@span}{c}{) suffix(})  					///
				span erepeat(\cmidrule(lr){@span}))  						///
			cells("mean(fmt(2) label(Mean)) sd(fmt(2) label(SD)) count(fmt(%9.0fc) label(Count)) min(fmt(%9.0fc) label(Min)) p50(fmt(%9.0fc) label(Median))  max(fmt(%9.0fc) label(Max))" ) ///
			label  nonumber noobs booktabs

		restore

		
*--------------
*		FIG A2
*--------------
		preserve
	
			est clear
			local X terms age age2 i.party_vv house0 

			forvalues i = 1(1)5 {
				label var sons_any_draft_running_b`i' "`i' year lead"

						
			*	1. 
			reghdfe vote_pro_cons  	sons_any_draft_running_b`i' 		///
				if kids_any_draft_running_b`i' == 1 , a(kids_num_now  sons_num_now) cluster(id_vv voteFE)			

						est sto reg1`i'		
						
			*	2. 
			reghdfe vote_pro_cons  	sons_any_draft_running_b`i' 	`X'	///
				if kids_any_draft_running_b`i' == 1 , a(kids_num_now  sons_num_now) cluster(id_vv voteFE)			
				
						est sto reg2`i'		

			*	3. 
			reghdfe vote_pro_cons  	sons_any_draft_running_b`i' 		///
				if kids_any_draft_running_b`i' == 1 , a(kids_num_now  sons_num_now voteFE ) cluster(id_vv voteFE)			

						est sto reg3`i'		
						
			*	4. 
			reghdfe vote_pro_cons  	sons_any_draft_running_b`i'  	`X'	/// 
				if kids_any_draft_running_b`i' == 1 , a(kids_num_now  sons_num_now voteFE ) cluster(id_vv voteFE)			
				
						est sto reg4`i'	
						
			}

			set scheme lean1

			coefplot (reg11) (reg12) (reg13) (reg14) (reg15) , ///
			keep(sons_any_draft_running_b*) xline(0) nokey offset(0) ///
			xtitle(Effect of Draft Age Son on Pro Draft Vote) name(a, replace) ci(90) title("(1) Baseline") msymbol(O) xlabel(-0.2 (0.05) 0.1)
			graph save "$output/JPE_REP_FA2a", replace  
			graph export "$output/JPE_REP_FA2a.png", as(png)   replace
			graph export "$output/JPE_REP_FA2a.pdf", as(pdf)   replace

			coefplot (reg21) (reg22) (reg23) (reg24) (reg25) , ///
			keep(sons_any_draft_running_b*) xline(0) nokey offset(0) ///
			xtitle(Effect of Draft Age Son on Pro Draft Vote) name(b, replace) ci(90) title("(2) Baseline + Controls") msymbol(O) xlabel(-0.2 (0.05) 0.1)
			graph save "$output/JPE_REP_FA2b", replace  
			graph export "$output/JPE_REP_FA2b.png", as(png)   replace
			graph export "$output/JPE_REP_FA2b.pdf", as(pdf)   replace

			coefplot  (reg31) (reg32) (reg33) (reg34) (reg35) , ///
			keep(sons_any_draft_running_b*) xline(0) nokey offset(0) ///
			xtitle(Effect of Draft Age Son on Pro Draft Vote) name(c, replace) ci(90) title("(3) Baseline + VFE") msymbol(O) xlabel(-0.2 (0.05) 0.1)
			graph save "$output/JPE_REP_FA2c", replace  
			graph export "$output/JPE_REP_FA2c.png", as(png)   replace
			graph export "$output/JPE_REP_FA2c.pdf", as(pdf)   replace

			coefplot  (reg41) (reg42) (reg43) (reg44) (reg45) , ///
			keep(sons_any_draft_running_b*) xline(0) nokey offset(0) ///
			xtitle(Effect of Draft Age Son on Pro Draft Vote) name(d, replace) ci(90) title("(4) Baseline + VFE + Controls") msymbol(O) xlabel(-0.2 (0.05) 0.1)
			graph save "$output/JPE_REP_FA2d", replace 
			graph export "$output/JPE_REP_FA2d.png", as(png)   replace
			graph export "$output/JPE_REP_FA2d.pdf", as(pdf)   replace
			
		restore
	
*--------------
*		FIG A3
*--------------
		preserve
		
			est clear
			local X terms age age2 i.party_vv house0 

			forvalues i = 1(1)5 {
				label var sons_any_draft_window_b`i' "`i' year lead"
						
			*	1. 
			reghdfe vote_pro_cons  	sons_any_draft_window_b`i' 		///
				if kids_any_draft_window_b`i' == 1 , a(kids_num_now  sons_num_now) cluster(id_vv voteFE)			

						est sto reg1`i'		
						
			*	2. 
			reghdfe vote_pro_cons  	sons_any_draft_window_b`i' 	`X'	///
				if kids_any_draft_window_b`i' == 1 , a(kids_num_now  sons_num_now) cluster(id_vv voteFE)			
				
						est sto reg2`i'		

			*	3. 
			reghdfe vote_pro_cons  	sons_any_draft_window_b`i' 		///
				if kids_any_draft_window_b`i' == 1 , a(kids_num_now  sons_num_now voteFE ) cluster(id_vv voteFE)			

						est sto reg3`i'		
						
			*	4. 
			reghdfe vote_pro_cons  	sons_any_draft_window_b`i'  	`X'	/// 
				if kids_any_draft_window_b`i' == 1 , a(kids_num_now  sons_num_now voteFE ) cluster(id_vv voteFE)			
				
						est sto reg4`i'	
						
			}

			set scheme lean1

			coefplot (reg11) (reg12) (reg13) (reg14) (reg15) , ///
			keep(sons_any_draft_window_b*) xline(0) nokey offset(0) ///
			xtitle(Effect of Draft Age Son on Pro Draft Vote) name(a, replace) ci(90) title("(1) Baseline") msymbol(O) xlabel(-0.2 (0.05) 0.1)
			graph save "$output/JPE_REP_FA3a", replace  
			graph export "$output/JPE_REP_FA3a.png", as(png)   replace
			graph export "$output/JPE_REP_FA3a.pdf", as(pdf)   replace

			coefplot (reg21) (reg22) (reg23) (reg24) (reg25) , ///
			keep(sons_any_draft_window_b*) xline(0) nokey offset(0) ///
			xtitle(Effect of Draft Age Son on Pro Draft Vote) name(b, replace) ci(90) title("(2) Baseline + Controls") msymbol(O) xlabel(-0.2 (0.05) 0.1)
			graph save "$output/JPE_REP_FA3b", replace  
			graph export "$output/JPE_REP_FA3b.png", as(png)   replace
			graph export "$output/JPE_REP_FA3b.pdf", as(pdf)   replace

			coefplot  (reg31) (reg32) (reg33) (reg34) (reg35) , ///
			keep(sons_any_draft_window_b*) xline(0) nokey offset(0) ///
			xtitle(Effect of Draft Age Son on Pro Draft Vote) name(c, replace) ci(90) title("(3) Baseline + VFE") msymbol(O) xlabel(-0.2 (0.05) 0.1)
			graph save "$output/JPE_REP_FA3c", replace  
			graph export "$output/JPE_REP_FA3c.png", as(png)   replace
			graph export "$output/JPE_REP_FA3c.pdf", as(pdf)   replace

			coefplot  (reg41) (reg42) (reg43) (reg44) (reg45) , ///
			keep(sons_any_draft_window_b*) xline(0) nokey offset(0) ///
			xtitle(Effect of Draft Age Son on Pro Draft Vote) name(d, replace) ci(90) title("(4) Baseline + VFE + Controls") msymbol(O) xlabel(-0.2 (0.05) 0.1)
			graph save "$output/JPE_REP_FA3d", replace 
			graph export "$output/JPE_REP_FA3d.png", as(png)   replace
			graph export "$output/JPE_REP_FA3d.pdf", as(pdf)   replace
		
		restore
		
	
*--------------
*		FIG A4 -- HAWKS AND DOVES
*--------------

	
*--------------
*		FIG A5 -- HAWKS AND DOVES
*--------------


*--------------
*		TABLE A5 
*--------------
		preserve
		
				est clear
				local X terms age age2 i.party_vv house0 
				
				label var sons_any_draft_running_lead4 	"\hspace{15pt} Son"
				label var sons_any_draft_running_b0 	"\hspace{15pt} Son"
				label var kids_any_draft_running_lead4 	"\hspace{15pt} Child"
				label var kids_any_draft_running_b0 	"\hspace{15pt} Child"
						
				
				
				reghdfe vote_pro_cons  	sons_any_draft_running_lead4 	///
										sons_any_draft_running_b0 		///
										kids_any_draft_running_lead4 	///
					if kids_any_draft_running_b4 == 1, a(kids_num_now  sons_num_now) cluster(id_vv voteFE)
								
							lincom sons_any_draft_running_b0 - kids_any_draft_running_lead4
							estadd sca est_sdkd = r(estimate)
							estadd sca se_sdkd = r(se)
							test sons_any_draft_running_b0 - kids_any_draft_running_lead4 = 0
							estadd local p_sdkd = trim("[`:display %9.2f r(p)']")
							
							test sons_any_draft_running_lead4 + kids_any_draft_running_lead4 = sons_any_draft_running_b0
							estadd local p_test = trim("[`:display %9.2f r(p)']")
							
					estadd local tst ""
					estadd local tst1 ""
					estadd local lfe "No"
					estadd local vfe "No"
					estadd local sfe "Yes"
					estadd local kfe "Yes"
					estadd local con "No"
					estadd sca clust1 = e(N_clust1)
					estadd sca clust2 = e(N_clust2)
					estadd ysumm
					
					est sto reg11
					
				reghdfe vote_pro_cons  	sons_any_draft_running_lead4 		///
										sons_any_draft_running_b0 			///
										kids_any_draft_running_lead4  `X' 	///
					if kids_any_draft_running_b4 == 1, a(kids_num_now  sons_num_now  ) cluster(id_vv voteFE)
								
							lincom sons_any_draft_running_b0 - kids_any_draft_running_lead4
							estadd sca est_sdkd = r(estimate)
							estadd sca se_sdkd = r(se)
							test sons_any_draft_running_b0 - kids_any_draft_running_lead4 = 0
							estadd local p_sdkd = trim("[`:display %9.2f r(p)']")
							
							test sons_any_draft_running_lead4 + kids_any_draft_running_lead4 = sons_any_draft_running_b0
							estadd local p_test = trim("[`:display %9.2f r(p)']")

					estadd local tst ""
					estadd local tst1 ""
					estadd local lfe "No"
					estadd local vfe "No"
					estadd local sfe "Yes"
					estadd local kfe "Yes"
					estadd local con "Yes"
					estadd sca clust1 = e(N_clust1)
					estadd sca clust2 = e(N_clust2)
					estadd ysumm
					
					est sto reg12
					
					
				reghdfe vote_pro_cons  	sons_any_draft_running_lead4 ///
										sons_any_draft_running_b0 ///
										kids_any_draft_running_lead4 ///
					if kids_any_draft_running_b4 == 1, a(kids_num_now  sons_num_now  voteFE) cluster(id_vv voteFE)
									
							lincom sons_any_draft_running_b0 - kids_any_draft_running_lead4
							estadd sca est_sdkd = r(estimate)
							estadd sca se_sdkd = r(se)
							test sons_any_draft_running_b0 - kids_any_draft_running_lead4 = 0
							estadd local p_sdkd = trim("[`:display %9.2f r(p)']")
							
							test sons_any_draft_running_lead4 + kids_any_draft_running_lead4 = sons_any_draft_running_b0
							estadd local p_test = trim("[`:display %9.2f r(p)']")
							
					estadd local tst ""
					estadd local tst1 ""
					estadd local lfe "No"
					estadd local vfe "Yes"
					estadd local sfe "Yes"
					estadd local kfe "Yes"
					estadd local con "No"
					estadd sca clust1 = e(N_clust1)
					estadd sca clust2 = e(N_clust2)
					estadd ysumm
					
					est sto reg13
					
				reghdfe vote_pro_cons  	sons_any_draft_running_lead4 ///
										sons_any_draft_running_b0 ///
										kids_any_draft_running_lead4 `X' ///
					if kids_any_draft_running_b4 == 1, a(kids_num_now  sons_num_now  voteFE) cluster(id_vv voteFE)
									
							lincom sons_any_draft_running_b0 - kids_any_draft_running_lead4
							estadd sca est_sdkd = r(estimate)
							estadd sca se_sdkd = r(se)
							test sons_any_draft_running_b0 - kids_any_draft_running_lead4 = 0
							estadd local p_sdkd = trim("[`:display %9.2f r(p)']")
							
							test sons_any_draft_running_lead4 + kids_any_draft_running_lead4 = sons_any_draft_running_b0
							estadd local p_test = trim("[`:display %9.2f r(p)']")
							
					estadd local tst ""
					estadd local tst1 ""
					estadd local lfe "No"
					estadd local vfe "Yes"
					estadd local sfe "Yes"
					estadd local kfe "Yes"
					estadd local con "Yes"
					estadd sca clust1 = e(N_clust1)
					estadd sca clust2 = e(N_clust2)
					estadd ysumm
					
					est sto reg14		
				
				
				esttab  reg11 reg12 reg13 reg14 using "$output/JPE_REP_TA5.tex", replace se noconstant drop(_cons* term* age* *party* house*) ///
				order("\underline{Leading Cohort} \vspace{-0.4cm}" sons_any_draft_running_lead4 kids_any_draft_running_lead4  "\underline{Current Cohort} \vspace{-0.4cm}" sons_any_draft_running_b0 ) ///
				star(* 0.100 ** 0.050 *** 0.010) b(%10.4f %10.4f)  nonotes   ///
				scalars(	"tst \underline{\emph{Additional Calculations}}" ///
							"est_sdkd \\ Current Son - Leading Child" "p_sdkd \hspace{15pt} p-value" ///
							"tst1 \\ Test: [Current Son - Leading Child] = [Leading Son]" ///
							"p_test \hspace{15pt}  p-value"  ///
							"con \hline \\ Controls" "vfe Vote FE"   ///
							"sfe Number of Sons FE" "kfe Number of Children FE"  /// 
							"clust1 Legislators" "clust2 Votes" ///
							"ymean Mean Dep. Var." "N Observations")  ///
				mgroups("Pro Draft Vote" , pattern(1 0 0 0)  ///
				prefix(\multicolumn{@span}{c}{) suffix(})  ///
				span erepeat(\cmidrule(lr){@span}))  ///
				sfmt( %~12s  %10.4f  %10.2f  %~12s   %10.2f %~12s %~12s  %~12s %~12s %~12s  %~12s  %10.2f %10.0f  )  /// 
				label obslast nomtitles booktabs		
			
	restore
	
	
*--------------
*		FIG A6 -- HAWKS AND DOVES
*--------------

	
*--------------
*		FIG A7 
*--------------
	rdrobust dem1 	son_RD  	if   	!missing(vote_pro_) &  sons_num_now ==1,  vce(cluster id_vv)  masspoints(adjust)
	rdplot 	 dem1 	son_RD 	if  	!missing(vote_pro_) &  sons_num_now ==1 & 	inrange(son_RD, - 6,  6), 			///
	h( 6) bw(mserd)  p(1)  c(0)  vce(cluster id_vv) graph_options(title("") ylabel(0 0.2 0.4 0.6 0.8 1 ) xlabel(-6(1)6)	 	///
	ytitle(Democrat) xtitle(Son's Age {&minus} Upper Cutoff) legend(off) ) masspoints(adjust)  ci(95) shade
	graph save "$output/JPE_REP_FA7a", replace 
	graph export "$output/JPE_REP_FA7a.png", as(png)   replace
	graph export "$output/JPE_REP_FA7a.pdf", as(pdf)   replace
		
	rdrobust age 		son_RD  	if   	!missing(vote_pro_) &  sons_num_now ==1 ,  vce(cluster id_vv)  masspoints(adjust)
	rdplot 	 age 		son_RD 	if  	!missing(vote_pro_) &  sons_num_now ==1 & 	inrange(son_RD, - 6,  6), 		///
	h( 6) bw(mserd)  p(1)  c(0)  vce(cluster id_vv) graph_options(title("") ylabel(50(5)70 ) xlabel(-6(1)6)	 				///
	ytitle(Age) xtitle(Son's Age {&minus} Upper Cutoff) legend(off) ) masspoints(adjust)  ci(95) shade
	graph save "$output/JPE_REP_FA7b", replace 
	graph export "$output/JPE_REP_FA7b.png", as(png)   replace
	graph export "$output/JPE_REP_FA7b.pdf", as(pdf)   replace
			
	rdrobust house0 		son_RD  	if   	!missing(vote_pro_) &  sons_num_now ==1,  vce(cluster id_vv)  masspoints(adjust)
	rdplot 	 house0 		son_RD 	if  	!missing(vote_pro_) &  sons_num_now ==1 & 	inrange(son_RD, - 6,  6), 	///
	h( 6) bw(mserd)  p(1)  c(0)  vce(cluster id_vv) graph_options(title("") ylabel(0 0.2 0.4 0.6 0.8 1 ) xlabel(-6(1)6)	 	///
	ytitle(Senator) xtitle(Son's Age {&minus} Upper Cutoff) legend(off) ) masspoints(adjust)  ci(95) shade
	graph save "$output/JPE_REP_FA7c", replace 
	graph export "$output/JPE_REP_FA7c.png", as(png)   replace
	graph export "$output/JPE_REP_FA7c.pdf", as(pdf)   replace
		
	rdrobust income 		son_RD  if   !missing(vote_pro_) &  sons_num_now ==1,  vce(cluster id_vv)  masspoints(adjust)
	rdplot 	 income 		son_RD 	if  	!missing(vote_pro_) &  sons_num_now ==1 & 	inrange(son_RD, - 6,  6), 	///
	h( 6) bw(mserd)  p(1)  c(0)  vce(cluster id_vv) graph_options(title("")  xlabel(-6(1)6)									///
	ytitle(Income) xtitle(Son's Age {&minus} Upper Cutoff) legend(off) ) masspoints(adjust)  ci(95) shade
	graph save "$output/JPE_REP_FA7d", replace 
	graph export "$output/JPE_REP_FA7d.png", as(png)   replace
	graph export "$output/JPE_REP_FA7d.pdf", as(pdf)   replace	
		
		
*--------------
*		FIG A8 -- HAWKS AND DOVES
*--------------	
	
	
*--------------
*		FIG A9 -- REPRODUCED FROM FORDHAM (2016)
*--------------	
		
	
*--------------
*		FIG A10
*--------------	
		est clear
		local X terms age age2 i.party_vv house0 

		foreach i in 1 2 5 {

		
			*	1
		reghdfe vote_pro_cons  	sons_any_draft_running_b4 ///
			if warFE`i' == 1 & kids_any_draft_running_b4 == 1, a(kids_num_now sons_num_now) cluster(id_vv voteFE)
						
			estadd local lfe "No"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg`i'1
 
 
 			*	2
		reghdfe vote_pro_cons  	sons_any_draft_running_b4 `X' ///
			if warFE`i' == 1 & kids_any_draft_running_b4 == 1, a(kids_num_now sons_num_now) cluster(id_vv voteFE)
 		
			estadd local lfe "No"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg`i'2
			
			
 			*	3
		reghdfe vote_pro_cons  	sons_any_draft_running_b4 ///
			if warFE`i' == 1 & kids_any_draft_running_b4 == 1, a(kids_num_now sons_num_now voteFE) cluster(id_vv voteFE)
						
			estadd local lfe "No"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg`i'3
 
 
 			*	4
		reghdfe vote_pro_cons  	sons_any_draft_running_b4 `X' ///
			if warFE`i' == 1 & kids_any_draft_running_b4 == 1, a(kids_num_now sons_num_now voteFE) cluster(id_vv voteFE)
 		
			estadd local lfe "No"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg`i'4
									
	}					
	
		set scheme lean1
		coefplot  	reg11 reg12 reg13 reg14, keep(sons_any_draft_running_b4) title("Impact of Draft Son on Pro Draft Votes by War") key  bylabel("World War I") ||		///
					reg21 reg22 reg23 reg24, keep(sons_any_draft_running_b4) bylabel("World War II") ||		///
					reg51 reg52 reg53 reg54, keep(sons_any_draft_running_b4) bylabel("Cold War") ||			///
					, byopts(compact cols(1)) xline(0) ///
		subtitle(, size(medium) margin(small) justification(left) color(black) bcolor(white) bmargin(top_bottom)) legend(rows(1))  msymbol(o) levels(90)

		/*
		graph save  "$output/JPE_REP_FA10.gph", replace
		graph export "$output/JPE_REP_FA10.png", as(png) replace
		graph export "$output/JPE_REP_FA10.pdf", as(pdf) replace
		*/
	
	
*--------------
*		TABLE A6
*--------------		
		*	A.
		est clear

		foreach i in dem1 house0 age terms income { 		
			reghdfe `i'  	sons_any_draft_running_b4 ///
				if kids_any_draft_running_b4 == 1, a(kids_num_now  sons_num_now) cluster(id_vv voteFE)
								
					estadd local sfe "Yes"
					estadd local kfe "Yes"
					estadd sca clust1 = e(N_clust1)
					estadd sca clust2 = e(N_clust2)
					estadd ysumm
					
					est sto regA`i'
		}
		esttab regA*  using  "$output/JPE_REP_TA6a.tex", replace se  drop(_cons) ///
		order( "\underline{\emph{Conditional Balance}} \vspace{-0.4cm}" sons_any_draft_running_b*) ///
		star(* 0.100 ** 0.050 *** 0.010) b(%10.3f %10.3f)  nonotes   ///
			scalars(	"sfe \\ Number of Sons FE" ///
						"kfe Number of Children FE"  /// 
						"clust1 Legislators" ///
						"clust2 Votes" ///
						"ymean Mean Dep. Var." ///
						"N Observations")  ///
			 sfmt(%~12s %~12s %~12s %~12s %10.3f %10.0f)  /// 
			 label obslast  booktabs			
					
					
		*	B.
		est clear

		foreach i in dem1 house0 age terms income { 
			reghdfe `i'  	sons_any_draft_running_b4 ///
				if kids_any_draft_running_b4 == 1, noa cluster(id_vv voteFE)

					estadd local sfe "No"
					estadd local kfe "No"
					estadd sca clust1 = e(N_clust1)
					estadd sca clust2 = e(N_clust2)
					estadd ysumm
					
					est sto regB`i'
		}
		esttab regB*  using  "$output/JPE_REP_TA6b.tex", replace se   ///
		order( "\underline{\emph{Unconditional Balance}} \vspace{-0.4cm}" sons_any_draft_running_b*) ///
		star(* 0.100 ** 0.050 *** 0.010) b(%10.3f %10.3f)  nonotes   ///
			scalars(	"sfe \\ Number of Sons FE" ///
						"kfe Number of Children FE"  /// 
						"clust1 Legislators" ///
						"clust2 Votes" ///
						"ymean Mean Dep. Var." ///
						"N Observations")  ///
			 sfmt(%~12s %~12s %~12s %~12s %10.3f %10.0f)  /// 
			 label obslast  booktabs		
					
			
*--------------
*		TABLE A7
*--------------	
		est clear
		local X terms age age2 i.party_vv house0 
		
		*	1.	
		reghdfe vote_pro_cons  	sons_any_draft_running_b4 ///
			if kids_any_draft_running_b4 == 1, noa cluster(id_vv voteFE)
						
			estadd local lfe "No"
			estadd local vfe "No"
			estadd local sfe "No"
			estadd local kfe "No"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg1
 
		*	2.	
		reghdfe vote_pro_cons  	sons_any_draft_running_b4 `X' ///
			if kids_any_draft_running_b4 == 1, noa cluster(id_vv voteFE)
 		
			estadd local lfe "No"
			estadd local vfe "No"
			estadd local sfe "No"
			estadd local kfe "No"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg2
			
			
		*	3.	
		reghdfe vote_pro_cons  	sons_any_draft_running_b4 ///
			if kids_any_draft_running_b4 == 1, a(voteFE) cluster(id_vv voteFE)
						
			estadd local lfe "No"
			estadd local vfe "Yes"
			estadd local sfe "No"
			estadd local kfe "No"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg3
 
		*	4.	
		reghdfe vote_pro_cons  	sons_any_draft_running_b4 `X' ///
			if kids_any_draft_running_b4 == 1, a(voteFE) cluster(id_vv voteFE)
 		
			estadd local lfe "No"
			estadd local vfe "Yes"
			estadd local sfe "No"
			estadd local kfe "No"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg4
			
					
		esttab reg1 reg2 reg3 reg4 using "$output/JPE_REP_TA7.tex", replace se noconstant keep(sons_any_draft_running_b*  ) order(sons_any_draft_running_b*  ) ///
			 star(* 0.100 ** 0.050 *** 0.010) b(%10.4f %10.4f)  nonotes   ///
			   scalars(	"con \\ Controls"			///
						"vfe Vote FE"  				/// 
						"sfe Number of Sons FE" 	///
						"kfe Number of Children FE" /// 
						"clust1 Legislators" 		///
						"clust2 Votes" 				///
						"ymean Mean Dep. Var." 		///
						"N Observations")  			///
				mgroups("Pro Draft Vote", pattern(1 0 0 0)  ///
				prefix(\multicolumn{@span}{c}{) suffix(})  ///
				span erepeat(\cmidrule(lr){@span}))  ///
			 sfmt(%~12s %~12s  %~12s %~12s %~12s  %~12s %10.2f %10.0f)  /// 
			 label obslast nomtitles booktabs
					
			
*--------------
*		TABLE A8
*--------------	
		preserve
			
			label var sons_any_draft_window_b2	"Draft Age Son"
			label var kids_any_draft_window_b2 	"Draft Age Child"
			label var sons_any_draft_rw_b42 	"Draft Age Son"
			label var kids_any_draft_rw_b42 	"Draft Age Child"
			label var sons_any_draft_rw_b4 		"Draft Age Son"
			label var kids_any_draft_rw_b4 		"Draft Age Child"

			est clear		
			*	1.
			reghdfe vote_pro_cons  	sons_any_draft_window_b2 ///
				if kids_any_draft_window_b2 == 1, a(kids_num_now  sons_num_now ) cluster(id_vv voteFE)
							
				estadd local lfe "No"
				estadd local vfe "No"
				estadd local sfe "Yes"
				estadd local kfe "Yes"
				estadd local con "No"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto reg1
	 
	 
			*	2.
			reghdfe vote_pro_cons  	sons_any_draft_window_b2 `X' ///
				if kids_any_draft_window_b2 == 1, a(kids_num_now  sons_num_now ) cluster(id_vv voteFE)
			
				estadd local lfe "No"
				estadd local vfe "No"
				estadd local sfe "Yes"
				estadd local kfe "Yes"
				estadd local con "Yes"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto reg2
				
				
			*	3.
			reghdfe vote_pro_cons  	sons_any_draft_window_b2 ///
				if kids_any_draft_window_b2 == 1, a(kids_num_now  sons_num_now voteFE) cluster(id_vv voteFE)
							
				estadd local lfe "No"
				estadd local vfe "Yes"
				estadd local sfe "Yes"
				estadd local kfe "Yes"
				estadd local con "No"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto reg3
	 
	 
			*	4.
			reghdfe vote_pro_cons  	sons_any_draft_window_b2 `X' ///
				if kids_any_draft_window_b2 == 1, a(kids_num_now  sons_num_now voteFE ) cluster(id_vv voteFE)
			
				estadd local lfe "No"
				estadd local vfe "Yes"
				estadd local sfe "Yes"
				estadd local kfe "Yes"
				estadd local con "Yes"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto reg4
				
			esttab  reg1 reg2 reg3 reg4 using "$output/JPE_REP_TA8a.tex", replace se noconstant keep(sons_any_draft_window_b2*  ) ///
				star(* 0.100 ** 0.050 *** 0.010) b(%10.4f %10.4f)  nonotes   ///
				scalars(	"con \\ Controls" ///
							"vfe Vote FE"   ///
							"sfe Number of Sons FE" ///
							"kfe Number of Children FE"  /// 
							"clust1 Legislators" ///
							"clust2 Votes" ///
							"ymean Mean Dep. Var." ///
							"N Observations")  ///
					mgroups("Pro Draft Vote: Window Votes (2Y Lead)" , pattern(1 0 0 0)  ///
					prefix(\multicolumn{@span}{c}{) suffix(})  ///
					span erepeat(\cmidrule(lr){@span}))  ///
				 sfmt(%~12s %~12s  %~12s %~12s %~12s  %~12s %10.2f %10.0f  )  /// 
				 label obslast nomtitles booktabs
						

			est clear		
			*	1.
			reghdfe vote_pro_cons  	sons_any_draft_rw_b42 ///
				if kids_any_draft_rw_b42 == 1, a(kids_num_now  sons_num_now ) cluster(id_vv voteFE)
							
				estadd local lfe "No"
				estadd local vfe "No"
				estadd local sfe "Yes"
				estadd local kfe "Yes"
				estadd local con "No"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto reg1
	 
	 
			*	2.
			reghdfe vote_pro_cons  	sons_any_draft_rw_b42 `X' ///
				if kids_any_draft_rw_b42 == 1, a(kids_num_now  sons_num_now ) cluster(id_vv voteFE)
			
				estadd local lfe "No"
				estadd local vfe "No"
				estadd local sfe "Yes"
				estadd local kfe "Yes"
				estadd local con "Yes"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto reg2
				
				
			*	3.
			reghdfe vote_pro_cons  	sons_any_draft_rw_b42 ///
				if kids_any_draft_rw_b42 == 1, a(kids_num_now  sons_num_now voteFE) cluster(id_vv voteFE)
							
				estadd local lfe "No"
				estadd local vfe "Yes"
				estadd local sfe "Yes"
				estadd local kfe "Yes"
				estadd local con "No"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto reg3
	 
	 
			*	4.
			reghdfe vote_pro_cons  	sons_any_draft_rw_b42 `X' ///
				if kids_any_draft_rw_b42 == 1, a(kids_num_now  sons_num_now voteFE ) cluster(id_vv voteFE)
			
				estadd local lfe "No"
				estadd local vfe "Yes"
				estadd local sfe "Yes"
				estadd local kfe "Yes"
				estadd local con "Yes"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto reg4
				
			esttab  reg1 reg2 reg3 reg4 using "$output/JPE_REP_TA8b.tex", replace se noconstant keep(sons_any_draft*  ) ///
				star(* 0.100 ** 0.050 *** 0.010) b(%10.4f %10.4f)  nonotes   ///
				scalars(	"con \\ Controls" ///
							"vfe Vote FE"   ///
							"sfe Number of Sons FE" ///
							"kfe Number of Children FE"  /// 
							"clust1 Legislators" ///
							"clust2 Votes" ///
							"ymean Mean Dep. Var." ///
							"N Observations")  ///
					mgroups("Pro Draft Vote: Window (2Y Lead) and Main (4Y Lead)" , pattern(1 0 0 0)  ///
					prefix(\multicolumn{@span}{c}{) suffix(})  ///
					span erepeat(\cmidrule(lr){@span}))  ///
				 sfmt(%~12s %~12s  %~12s %~12s %~12s  %~12s %10.2f %10.0f  )  /// 
				 label obslast nomtitles booktabs
				 

			est clear		
			*	1.
			reghdfe vote_pro_cons  	sons_any_draft_rw_b4 ///
				if kids_any_draft_rw_b4 == 1, a(kids_num_now  sons_num_now ) cluster(id_vv voteFE)
							
				estadd local lfe "No"
				estadd local vfe "No"
				estadd local sfe "Yes"
				estadd local kfe "Yes"
				estadd local con "No"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto reg1
	 
	 
			*	2.
			reghdfe vote_pro_cons  	sons_any_draft_rw_b4 `X' ///
				if kids_any_draft_rw_b4 == 1, a(kids_num_now  sons_num_now ) cluster(id_vv voteFE)
			
				estadd local lfe "No"
				estadd local vfe "No"
				estadd local sfe "Yes"
				estadd local kfe "Yes"
				estadd local con "Yes"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto reg2
				
				
			*	3.
			reghdfe vote_pro_cons  	sons_any_draft_rw_b4 ///
				if kids_any_draft_rw_b4 == 1, a(kids_num_now  sons_num_now voteFE) cluster(id_vv voteFE)
							
				estadd local lfe "No"
				estadd local vfe "Yes"
				estadd local sfe "Yes"
				estadd local kfe "Yes"
				estadd local con "No"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto reg3
	 
	 
			*	4.
			reghdfe vote_pro_cons  	sons_any_draft_rw_b4 `X' ///
				if kids_any_draft_rw_b4 == 1, a(kids_num_now  sons_num_now voteFE ) cluster(id_vv voteFE)
			
				estadd local lfe "No"
				estadd local vfe "Yes"
				estadd local sfe "Yes"
				estadd local kfe "Yes"
				estadd local con "Yes"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto reg4
				
			esttab  reg1 reg2 reg3 reg4 using "$output/JPE_REP_TA8c.tex", replace se noconstant keep(sons_any_draft* ) ///
				star(* 0.100 ** 0.050 *** 0.010) b(%10.4f %10.4f)  nonotes   ///
				scalars(	"con \\ Controls" ///
							"vfe Vote FE"   ///
							"sfe Number of Sons FE" ///
							"kfe Number of Children FE"  /// 
							"clust1 Legislators" ///
							"clust2 Votes" ///
							"ymean Mean Dep. Var." ///
							"N Observations")  ///
					mgroups("Pro Draft Vote: Window (4Y Lead) and Main (4Y Lead)" , pattern(1 0 0 0)  ///
					prefix(\multicolumn{@span}{c}{) suffix(})  ///
					span erepeat(\cmidrule(lr){@span}))  ///
				 sfmt(%~12s %~12s  %~12s %~12s %~12s  %~12s %10.2f %10.0f  )  /// 
				 label obslast nomtitles booktabs
					
		restore
			
			
*--------------
*		TABLE A9 -- HAWKS AND DOVES
*--------------	
	
	
*--------------
*		TABLE A10
*--------------	
		est clear		
		*	1.		
		rdrobust vote_pro_ 	son_RD  if  sons_num_now ==1 , all  vce(cluster id_vv)  masspoints(adjust)
				estadd local rv "Son"
				estadd local fe "No"
				estadd local con "No"
				estadd sca clust1 = e(N_clust1)
				est sto regRD1
				
		*	2.		
		rdrobust vote_pro_ 	dtr_RD  if  dtrs_num_now ==1, all  vce(cluster id_vv)  masspoints(adjust)
				estadd local rv "Daughter"
				estadd local fe "No"
				estadd local con "No"
				estadd sca clust1 = e(N_clust1)
				est sto regRD2

		*	3.		
		rdrobust dem1 		son_RD  if   !missing(vote_pro_) &  sons_num_now ==1,  all vce(cluster id_vv)  masspoints(adjust)
				estadd local rv "Son"
				estadd local fe "No"
				estadd local con "No"
				estadd sca clust1 = e(N_clust1)
				est sto regRD3
				
		*	4.						
		rdrobust age 		son_RD  if   !missing(vote_pro_) &  sons_num_now ==1, all vce(cluster id_vv)  masspoints(adjust)
				estadd local rv "Son"
				estadd local fe "No"
				estadd local con "No"
				estadd sca clust1 = e(N_clust1)
				est sto regRD4

		*	5.						
		rdrobust house0 		son_RD  if   !missing(vote_pro_) &  sons_num_now ==1, all  vce(cluster id_vv)  masspoints(adjust)
				estadd local rv "Son"
				estadd local fe "No"
				estadd local con "No"
				estadd sca clust1 = e(N_clust1)
				est sto regRD5

		*	6.			
		rdrobust income 		son_RD  if   !missing(vote_pro_) &  sons_num_now ==1, all  vce(cluster id_vv)  masspoints(adjust)
				estadd local rv "Son"
				estadd local fe "No"
				estadd local con "No"
				estadd sca clust1 = e(N_clust1)
				est sto regRD6

		esttab regRD1 regRD2 regRD3 regRD4 regRD5 regRD6  using "$output/JPE_REP_TA10.tex", replace se ///
			star(* 0.100 ** 0.050 *** 0.010) b(%10.4f %10.4f)  nonotes   ///
			scalars(	"rv  \\ Running Variable: Age of" ///
						"h_l  Bandwidth"  )  ///
				mgroups("Pro Draft Vote" "Democrat" "Age" "Senator" "Income", pattern(1 0 1 1 1 1)  ///
				prefix(\multicolumn{@span}{c}{) suffix(})  ///
			span erepeat(\cmidrule(lr){@span}))  ///
			sfmt(  %~12s %10.2f   %10.0f  )  /// 
			label obslast nomtitles varlabels(Conventional "I(Age > Upper Cutoff)" ///
				Bias-corrected "\addlinespace \hspace{0.3cm} --- CCT Bias-Corrected Estimate" ///
				Robust "\addlinespace \hspace{0.3cm} --- with CCT Robust SE") booktabs
	 
	 
*--------------
*		TABLE A11
*--------------	
		est clear

		cap tab voteFE, gen(FE_vote)
		cap tab party_vv, gen (FE_party)

		*1
		rdrobust vote_pro_ 	son_RD  if  sons_num_now ==1 ,  vce(cluster id_vv)  masspoints(adjust) all  covs(terms age age2 house0 FE_party*)  covs_drop(pinv)
				estadd local rv "Son"
				estadd local fe "No"
				estadd local con "Yes"
				estadd sca clust1 = e(N_clust1)
				est sto regRD1
				
		*2
		rdrobust vote_pro_ 	dtr_RD  if  dtrs_num_now ==1,  vce(cluster id_vv)  masspoints(adjust) all  covs(terms age age2 house0 FE_party*) covs_drop(pinv)
				estadd local rv "Daughter"
				estadd local fe "No"
				estadd local con "Yes"
				estadd sca clust1 = e(N_clust1)
				est sto regRD2
				
		*3		
		rdrobust vote_pro_ 		son_RD  if     sons_num_now ==1,  vce(cluster id_vv)  masspoints(adjust) all  covs(FE_vote*) covs_drop(pinv)
				estadd local rv "Son"
				estadd local fe "Yes"
				estadd local con "No"
				estadd sca clust1 = e(N_clust1)
				est sto regRD3
				
		*4
		rdrobust vote_pro_ 		dtr_RD  if    sons_num_now ==1,  vce(cluster id_vv)  masspoints(adjust) all  covs(FE_vote*) covs_drop(pinv)
				estadd local rv "Daughter"
				estadd local fe "Yes"
				estadd local con "No"
				estadd sca clust1 = e(N_clust1)
				est sto regRD4
				
		*5			
		rdrobust vote_pro_ 		son_RD  if    sons_num_now ==1,  vce(cluster id_vv)  masspoints(adjust) all  covs(terms age age2  FE_party* FE_vote*) covs_drop(pinv)
				estadd local rv "Son"
				estadd local fe "Yes"
				estadd local con "Yes"
				estadd sca clust1 = e(N_clust1)
				est sto regRD5
				
		*6					
		rdrobust vote_pro_ 		dtr_RD  if     sons_num_now ==1,  vce(cluster id_vv)  masspoints(adjust) all  covs(terms age age2  FE_party* FE_vote*) covs_drop(pinv)
				estadd local rv "Daughter"
				estadd local fe "Yes"
				estadd local con "Yes"
				estadd sca clust1 = e(N_clust1)
				est sto regRD6

		esttab regRD1 regRD2 regRD3 regRD4 regRD5 regRD6  using "$output/JPE_REP_TA11.tex", replace se ///
			star(* 0.100 ** 0.050 *** 0.010) b(%10.4f %10.4f)  nonotes   ///
			scalars(	"rv  \\ Running Variable: Age of" ///
						"h_l  Bandwidth" ///
						"con  Controls" ///
						"fe  Vote FE"  )  ///
			mgroups("Pro Draft Vote", pattern(1 0 0 0 0 0 )  ///
			prefix(\multicolumn{@span}{c}{) suffix(})  ///
			span erepeat(\cmidrule(lr){@span}))  ///
			sfmt( %~12s %10.2f  %~12s %~12s  %10.0f  )  /// 
			label obslast nomtitles varlabels(Conventional "I(Age > Upper Cutoff)" ///
				Bias-corrected "\addlinespace \hspace{0.3cm} --- CCT Bias-Corrected Estimate" ///
				Robust "\addlinespace \hspace{0.3cm} --- with CCT Robust SE") booktabs
	 
	 
*--------------
*		TABLE A12 -- HAWKS AND DOVES
*--------------	


*--------------
*		TABLE A13
*--------------	
		est clear
		set graphics off
		
		*1
		qui eventdd	vote_pro_cons	i.youngest27_kid_event*,	timevar(youngest_son_event )	leads(5) lags(5) accum ///
		method(hdfe, absorb(id_vv   kids_num_now  sons_num_now) cluster(id_vv voteFE)) 
		 
					estadd local lfe "Yes"
					estadd local vfe "No"
					estadd local sfe "Yes"
					estadd local kfe "Yes"
					estadd local con "No"
					estadd sca clust1 = e(N_clust1)
					estadd sca clust2 = e(N_clust2)
					estadd ysumm

					est sto reg01		
		
		*2
		qui eventdd vote_pro_cons   i.youngest27_kid_event* terms age age2 i.party_vv i.house0, 	timevar(youngest_son_event ) leads(5) lags(5) accum ///
		method(hdfe, absorb(id_vv   kids_num_now  sons_num_now) cluster(id_vv voteFE)) 
	
					estadd local lfe "Yes"
					estadd local vfe "No"
					estadd local sfe "Yes"
					estadd local kfe "Yes"
					estadd local con "Yes"
					estadd sca clust1 = e(N_clust1)
					estadd sca clust2 = e(N_clust2)
					estadd ysumm

					est sto reg02
				
		*3	
		qui eventdd vote_pro_cons   i.youngest27_kid_event*,	timevar(youngest_son_event ) leads(5) lags(5) accum ///
		method(hdfe, absorb(id_vv voteFE  kids_num_now  sons_num_now) cluster(id_vv voteFE)) 
 

					estadd local lfe "Yes"
					estadd local vfe "Yes"
					estadd local sfe "Yes"
					estadd local kfe "Yes"
					estadd local con "No"
					estadd sca clust1 = e(N_clust1)
					estadd sca clust2 = e(N_clust2)
					estadd ysumm

					est sto reg03
				
		*4
		qui eventdd vote_pro_cons   i.youngest27_kid_event* age terms age2 i.party_vv i.house0, 	timevar(youngest_son_event ) leads(5) lags(5) accum ///
		method(hdfe, absorb(id_vv voteFE  kids_num_now  sons_num_now) cluster(id_vv voteFE))
	
					estadd local lfe "Yes"
					estadd local vfe "Yes"
					estadd local sfe "Yes"
					estadd local kfe "Yes"
					estadd local con "Yes"
					estadd sca clust1 = e(N_clust1)
					estadd sca clust2 = e(N_clust2)
					estadd ysumm

					est sto reg04
					
		set graphics on
		
		esttab reg01 reg02  reg03 reg04  using "$output/JPE_REP_TA13.tex", replace se noconstant keep(l* ) order(l* ) ///
			star(* 0.100 ** 0.050 *** 0.010) b(%10.4f %10.4f)  nonotes   ///
			scalars(	"con \\ Controls" ///
						"vfe Vote FE" ///
						"lfe Legislator FE"  ///
						"sfe Number of Sons FE" ///
						"kfe Number of Children FE"  /// 
						"clust1 Legislators" ///
						"clust2 Votes" ///
						"ymean Mean Dep. Var." ///
						"N Observations")  ///
		mgroups("Pro Draft Vote" , pattern(1 0 0 0)  ///
		prefix(\multicolumn{@span}{c}{) suffix(})  ///
		span erepeat(\cmidrule(lr){@span}))  ///
		sfmt(%~12s %~12s %~12s %~12s %~12s %~12s %~12s %10.2f %10.0f)  /// 
		label obslast nomtitles booktabs	///
		varlabels(	lead5 "Son Relative Age <= - 5" ///
					lead4 "Son Relative Age = - 4" ///
					lead3 "Son Relative Age = - 3" ///
					lead2 "Son Relative Age = - 2" ///
					lag0 "Son Relative Age = 0" ///
					lag1 "Son Relative Age = +1" ///
					lag2 "Son Relative Age = +2" ///
					lag3 "Son Relative Age = +3" ///
					lag4 "Son Relative Age = +4" ///
					lag5 "Son Relative Age >= +5")
				
				
*--------------
*		TABLE A14
*--------------	
		est clear
		local X terms age age2 i.party_vv house0 
		
		foreach i in 1 2 5 {
		
			*	1.		
			reghdfe vote_pro_cons  	electionyear `X'    ///
				if warFE`i' == 1 & !missing(sons_any_draft_running_b4), a(voteFE) cluster(id_vv voteFE)	
				
				estadd local lfe "No"
				estadd local vfe "Yes"
				estadd local sfe "No"
				estadd local kfe "No"
				estadd local con "Yes"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto reg`i'1	
		 
				
			reghdfe vote_pro_cons  	electionyear  `X'    ///
				if warFE`i' == 1 & !missing(sons_any_draft_running_b4), a(id_vv voteFE) cluster(id_vv voteFE)
									
				estadd local lfe "Yes"
				estadd local vfe "Yes"
				estadd local sfe "No"
				estadd local kfe "No"
				estadd local con "Yes"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto reg`i'2	
						
	}

		esttab  reg11 reg12  reg21 reg22  reg51 reg52 using "$output/JPE_REP_TA14_main.tex", replace se noconstant keep(electionyear) ///
			star(* 0.100 ** 0.050 *** 0.010) b(%10.4f %10.4f)  nonotes   ///
			scalars(	"lfe \\ Legislator FE" ///
						"vfe  Vote FE"  ///
						"sfe Number of Sons FE" ///
						"kfe Number of Children FE" ///
						"con Other Controls" /// 
						"clust1 Legislators"  ///
						"clust2 Votes" ///
						"ymean Mean Dep. Var." ///
						"N Observations")  ///
			mgroups("WW1" "WW2" "Cold War" , pattern(1  0 1 0 1 0)  ///
			prefix(\multicolumn{@span}{c}{) suffix(})  ///
			span erepeat(\cmidrule(lr){@span}))  ///
			sfmt(%~12s %~12s  %~12s %~12s %~12s  %~12s %~12s  ///
			%10.3f %10.0f  )  /// 
			label obslast nomtitles booktabs

						
*--------------
*		TABLE A15 -- HAWKS AND DOVES
*--------------	


*--------------
*		TABLE A16 -- HAWKS AND DOVES
*--------------	


*--------------
*		TABLE A17 -- ABOVE
*--------------	


*--------------
*		TABLE A18
*--------------	
		est clear		
		foreach i in 1 2 5 {
							
			reghdfe win  sons_any_draft_running_b4 	///
				if warFE`i' == 1 & kids_any_draft_running_b4 == 1, a(kids_num_now  sons_num_now) cluster(id_vv voteFE)	
	 		
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			est sto reg`i'1
					
			reghdfe win  sons_any_draft_running_b4 	///
				if warFE`i' == 1 & kids_any_draft_running_b4 == 1, a(kids_num_now  sons_num_now voteFE) cluster(id_vv voteFE)	
	 		
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg`i'2							
	}
	
		esttab  reg11 reg12  reg21 reg22  reg51 reg52  using "$output/JPE_REP_TA18.tex", replace se noconstant keep(sons_any_draft_running_b*) order(sons_any_draft_running_b* ) ///
		star(* 0.100 ** 0.050 *** 0.010) b(%10.4f %10.4f)  nonotes   ///
		scalars(	"vfe \\ Vote FE"  			///
					"sfe Number of Sons FE" 	///
					"kfe Number of Children FE" ///
					"con Other Controls" 		/// 
					"clust1 Legislators"  		///
					"clust2 Votes" 				///
					"ymean Mean Dep. Var." 		///
					"N Observations")  			///
		mgroups("WW1" "WW2" "Cold War" , pattern(1 0  1 0  1  0)  	///
		prefix(\multicolumn{@span}{c}{) suffix(})  					///
		span erepeat(\cmidrule(lr){@span}))  						///
		sfmt(%~12s  %~12s %~12s %~12s  %~12s  %~12s %10.3f %10.0f)  	/// 
		label obslast nomtitles booktabs


*--------------
*		TABLE A19
*--------------		 
		est clear
		local X terms age age2 i.party_vv house0 
		
			reghdfe vote_pro_cons  	sons_any_draft_running_b4 		///
									kids_any_draft_running_b4 `X' 	///
									presparty, a(kids_num_now  sons_num_now  ) cluster(id_vv voteFE)	
					
					estadd local lfe "No"
					estadd local vfe "No"
					estadd local sfe "Yes"
					estadd local kfe "Yes"
					estadd local con "Yes"
					estadd sca clust1 = e(N_clust1)
					estadd sca clust2 = e(N_clust2)
					estadd ysumm
					
					est sto reg11
			 
			reghdfe vote_pro_cons  	sons_any_draft_running_b4 		///
									kids_any_draft_running_b4 `X' 	///
									presparty, a(kids_num_now sons_num_now voteFE ) cluster(id_vv voteFE)	
					
					estadd local lfe "No"
					estadd local vfe "Yes"
					estadd local sfe "Yes"
					estadd local kfe "Yes"
					estadd local con "Yes"
					estadd sca clust1 = e(N_clust1)
					estadd sca clust2 = e(N_clust2)
					estadd ysumm
					
					est sto reg12	 
			 
			 
			reghdfe vote_pro_cons  	sons_any_draft_running_b4 		///
									kids_any_draft_running_b4 `X' 	///
									presparty, a(kids_num_now  sons_num_now id_vv  ) cluster(id_vv voteFE)	
					
					estadd local lfe "Yes"
					estadd local vfe "No"
					estadd local sfe "Yes"
					estadd local kfe "Yes"
					estadd local con "Yes"
					estadd sca clust1 = e(N_clust1)
					estadd sca clust2 = e(N_clust2)
					estadd ysumm
					
					est sto reg13	 
			 
			 
			reghdfe vote_pro_cons  	sons_any_draft_running_b4 		///
									kids_any_draft_running_b4 `X' 	///
									presparty, a(kids_num_now  sons_num_now id_vv voteFE ) cluster(id_vv voteFE)
									
					estadd local lfe "Yes"
					estadd local vfe "Yes"
					estadd local sfe "Yes"
					estadd local kfe "Yes"
					estadd local con "Yes"
					estadd sca clust1 = e(N_clust1)
					estadd sca clust2 = e(N_clust2)
					estadd ysumm
					
					est sto reg14						
								
		esttab reg11 reg12 reg13 reg14  using "$output/JPE_REP_TA19.tex", replace se noconstant keep(sons_any_draft_running_b* kids_any_draft_running_b* presparty) ///
		order(sons_any_draft_running_b* slead* kids_any_draft_running_b* presparty ) ///
		star(* 0.100 ** 0.050 *** 0.010) b(%10.4f %10.4f)  nonotes   ///
		scalars(	"lfe \\ Legislator FE" "vfe  Vote FE"  ///
					"sfe Number of Sons FE" ///
					"kfe Number of Children FE" ///
					"con Other Controls" /// 
					"clust1 Legislators" ///
					"clust2 Votes" ///
					"ymean Mean Dep. Var." ///
					"N Observations")  ///
		mgroups("Pro Draft Vote" , pattern(1 0 0 0)  ///
		prefix(\multicolumn{@span}{c}{) suffix(})  ///
		span erepeat(\cmidrule(lr){@span}))  ///
		sfmt(%~12s %~12s  %~12s %~12s %~12s  %~12s %~12s  ///
		%10.3f %10.0f  )  /// 
		label obslast nomtitles booktabs
			
			
			
***************************************************
***************************************************
***************************************************
*	B.	REPLICATION CODE: HAWKS AND DOVES
***************************************************
***************************************************
***************************************************

*	Load data
use "mhm_nokin_JPE_REP_DATA.dta", clear


*--------------
*		FIG 4a 
*--------------	

		est clear
		local X terms age age2 i.party_vv house0 
		
		foreach i in 1 2 5 {
		
			reghdfe vote_pro_cons  	electionyear `X'    ///
				if warFE`i' == 1 & !missing(sons_any_draft_running_b4), a(voteFE) cluster(id_vv voteFE)	
				
				estadd local lfe "No"
				estadd local vfe "Yes"
				estadd local sfe "No"
				estadd local kfe "No"
				estadd local con "Yes"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto regM`i'1	
		 
				
			reghdfe vote_pro_cons  	electionyear  `X'    ///
				if warFE`i' == 1 & !missing(sons_any_draft_running_b4), a(id_vv voteFE) cluster(id_vv voteFE)
									
				estadd local lfe "Yes"
				estadd local vfe "Yes"
				estadd local sfe "No"
				estadd local kfe "No"
				estadd local con "Yes"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto regM`i'2	
						
									
			reghdfe w_h1d0_mode  	electionyear `X' ///
				if warFE`i' == 1 , a(voteFE ) cluster(id_vv voteFE)	
					
					estadd local lfe "No"
					estadd local vfe "Yes"
					estadd local sfe "No"
					estadd local kfe "No"
					estadd local con "Yes"
					estadd sca clust1 = e(N_clust1)
					estadd sca clust2 = e(N_clust2)
					estadd ysumm
					
					est sto regS`i'1
			 
					
			reghdfe w_h1d0_mode  	electionyear `X' ///
				if warFE`i' == 1 , a(id_vv  voteFE ) cluster(id_vv voteFE)	
								
					estadd local lfe "Yes"
					estadd local vfe "Yes"
					estadd local sfe "No"
					estadd local kfe "No"
					estadd local con "Yes"
					estadd sca clust1 = e(N_clust1)
					estadd sca clust2 = e(N_clust2)
					estadd ysumm
					
					est sto regS`i'2
								
			}

		set scheme lean1	
		coefplot 	regM11 regS11  regM12 regS12  , keep(electionyear)   	bylabel("World War I")  || 		///
					regM21 regS21  regM22 regS22  , keep(electionyear)  	bylabel("World War II") || 		///
					regM51 regS51  regM52 regS52  , keep(electionyear)  	bylabel("Cold War") 	||		///
					,  byopts(compact cols(1))  ///
		subtitle(, size(medium) margin(small) justification(left) color(black) bcolor(white) bmargin(top_bottom)) legend(rows(1)) ///
		xtitle("Impact of Election Year on Pro Draft or Hawkish Vote")  levels(90) xline(0)

		/*
		graph save  "$output/JPE_REP_F4a.gph", replace
		graph export "$output/JPE_REP_F4a.eps", as(eps) preview(off) replace
		graph export "$output/JPE_REP_F4a.png", as(png) replace
		graph export "$output/JPE_REP_F4a.pdf", as(pdf) replace
		*/


*--------------
*		FIG A4 & A5
*--------------
		est clear
		
		preserve
		
			local X terms age age2 i.party_vv house0 
							
			foreach y of varlist   w_h1d0_mode w_h1d0_extrap_mode   {

				forvalues i = 1(1)5 {
					
					 label var sons_any_draft_ar_b`i' "`i' year lead"
		
				*			1.	KFE SFE
					reghdfe `y'	sons_any_draft_ar_b`i' ///
						if kids_any_draft_ar_b`i' == 1 &  all_draft == 1, a(kids_num_now  sons_num_now ) cluster(id_vv voteFE)
											
								estadd local lfe "No"
								estadd local vfe "No"
								estadd local sfe "Yes"
								estadd local kfe "Yes"
								estadd local con "No"
								estadd sca clust1 = e(N_clust1)
								estadd sca clust2 = e(N_clust2)
								estadd ysumm
								
								est sto reg1`i' 
					 
					reghdfe `y'	sons_any_draft_ar_b`i' `X' ///
						if kids_any_draft_ar_b`i' == 1 &  all_draft == 1, a(kids_num_now  sons_num_now ) cluster(id_vv voteFE)
							
								estadd local lfe "No"
								estadd local vfe "No"
								estadd local sfe "Yes"
								estadd local kfe "Yes"
								estadd local con "Yes"
								estadd sca clust1 = e(N_clust1)
								estadd sca clust2 = e(N_clust2)
								estadd ysumm
								
								est sto reg2`i'			
								
					*			2.	VFE KFE SFE
					reghdfe `y'   	sons_any_draft_ar_b`i' ///
						if kids_any_draft_ar_b`i' == 1 &  all_draft == 1, a(kids_num_now  sons_num_now voteFE) cluster(id_vv voteFE)
											
								estadd local lfe "No"
								estadd local vfe "Yes"
								estadd local sfe "Yes"
								estadd local kfe "Yes"
								estadd local con "No"
								estadd sca clust1 = e(N_clust1)
								estadd sca clust2 = e(N_clust2)
								estadd ysumm
								
								est sto reg3`i'			
								 
					reghdfe `y'   	sons_any_draft_ar_b`i' `X' ///
						if kids_any_draft_ar_b`i' == 1 &  all_draft == 1, a(kids_num_now  sons_num_now voteFE ) cluster(id_vv voteFE)
							
								estadd local lfe "No"
								estadd local vfe "Yes"
								estadd local sfe "Yes"
								estadd local kfe "Yes"
								estadd local con "Yes"
								estadd sca clust1 = e(N_clust1)
								estadd sca clust2 = e(N_clust2)
								estadd ysumm
								
								est sto reg4`i'				
		
	}
			
			set scheme lean1

			coefplot (reg11) (reg12) (reg13) (reg14) (reg15) , keep(sons_any_draft_ar_b1 sons_any_draft_ar_b2 sons_any_draft_ar_b3 sons_any_draft_ar_b4 sons_any_draft_ar_b5) ///
			xline(0) nokey offset(0) xtitle(Effect of Draft Age Son on Pro Draft Vote) name(c, replace) ci(90   ) title("(1) Baseline") msymbol(O) xlabel(-0.2 (0.05) 0.1)
			graph save "$output/JPE_REP_FA4A5_`y'_a", replace  // 102 clusters
			graph export "$output/JPE_REP_FA4A5_`y'_a.png", as(png)   replace
			graph export "$output/JPE_REP_FA4A5_`y'_a.pdf", as(pdf)   replace

			coefplot (reg21) (reg22) (reg23) (reg24) (reg25) , keep(sons_any_draft_ar_b1 sons_any_draft_ar_b2 sons_any_draft_ar_b3 sons_any_draft_ar_b4 sons_any_draft_ar_b5) ///
			xline(0) nokey offset(0) xtitle(Effect of Draft Age Son on Pro Draft Vote) name(d, replace) ci(90   ) title("(2) Baseline + Controls") msymbol(O) xlabel(-0.2 (0.05) 0.1)
			graph save "$output/JPE_REP_FA4A5_`y'_b", replace  // 102 clusters
			graph export "$output/JPE_REP_FA4A5_`y'_b.png", as(png)   replace
			graph export "$output/JPE_REP_FA4A5_`y'_b.pdf", as(pdf)   replace

			coefplot (reg31) (reg32) (reg33) (reg34) (reg35) , keep(sons_any_draft_ar_b1 sons_any_draft_ar_b2 sons_any_draft_ar_b3 sons_any_draft_ar_b4 sons_any_draft_ar_b5) ///
			xline(0) nokey offset(0) xtitle(Effect of Draft Age Son on Pro Draft Vote) name(e, replace) ci(90   ) title("(3) Baseline + VFE") msymbol(O) xlabel(-0.2 (0.05) 0.1)
			graph save "$output/JPE_REP_FA4A5_`y'_c", replace  // 102 clusters
			graph export "$output/JPE_REP_FA4A5_`y'_c.png", as(png)   replace
			graph export "$output/JPE_REP_FA4A5_`y'_c.pdf", as(pdf)   replace

			coefplot (reg41) (reg42) (reg43) (reg44) (reg45) , keep(sons_any_draft_ar_b1 sons_any_draft_ar_b2 sons_any_draft_ar_b3 sons_any_draft_ar_b4 sons_any_draft_ar_b5) ///
			xline(0) nokey offset(0) xtitle(Effect of Draft Age Son on Pro Draft Vote) name(f, replace) ci(90   ) title("(4) Baseline + Controls + VFE") msymbol(O) xlabel(-0.2 (0.05) 0.1)
			graph save "$output/JPE_REP_FA4A5_`y'_d", replace  // 102 clusters
			graph export "$output/JPE_REP_FA4A5_`y'_d.png", as(png)   replace
			graph export "$output/JPE_REP_FA4A5_`y'_d.pdf", as(pdf)   replace	

			}

		restore

	
*--------------
*		TABLE A9
*--------------	
	

*		A. All Draft Votes

		est clear
		local X terms age age2 i.party_vv house0 

			reghdfe w_h1d0_mode   	sons_any_draft_arw_b4 ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 1, a(kids_num_now  sons_num_now ) cluster(id_vv voteFE)
						
			estadd local lfe "No"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg111
 
 
			reghdfe w_h1d0_mode   	sons_any_draft_arw_b4 `X' ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 1, a(kids_num_now  sons_num_now ) cluster(id_vv voteFE)
 		
			estadd local lfe "No"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg112
			
			
			reghdfe w_h1d0_mode  	sons_any_draft_arw_b4 ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 1, a(kids_num_now  sons_num_now voteFE) cluster(id_vv voteFE)
						
			estadd local lfe "No"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg113
 
 
			reghdfe w_h1d0_mode  	sons_any_draft_arw_b4  `X' ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 1, a(kids_num_now  sons_num_now voteFE ) cluster(id_vv voteFE)
 		
			estadd local lfe "No"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg114
	

			reghdfe w_h1d0_extrap_mode   	sons_any_draft_arw_b4 ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 1, a(kids_num_now  sons_num_now ) cluster(id_vv voteFE)
						
			estadd local lfe "No"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg121
 
 
			reghdfe w_h1d0_extrap_mode   	sons_any_draft_arw_b4  `X' ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 1, a(kids_num_now  sons_num_now ) cluster(id_vv voteFE)
 		
			estadd local lfe "No"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg122
			
			
			reghdfe w_h1d0_extrap_mode  	sons_any_draft_arw_b4 ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 1, a(kids_num_now  sons_num_now voteFE) cluster(id_vv voteFE)
						
			estadd local lfe "No"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg123
 
 
			reghdfe w_h1d0_extrap_mode  	sons_any_draft_arw_b4  `X' ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 1, a(kids_num_now  sons_num_now voteFE ) cluster(id_vv voteFE)
 		
			estadd local lfe "No"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg124
	
			esttab   reg111 reg112 reg113 reg114 reg121 reg122 reg123 reg124 using "$output/JPE_REP_TA9a.tex", replace se  drop(terms age age2 terms age age2 *part* *cons* *house*  ) ///
			order("\underline{\emph{All Draft Votes}} \vspace{-0.4cm}" *) ///
			star(* 0.100 ** 0.050 *** 0.010) b(%10.4f %10.4f)  nonotes   ///
			scalars(	"con \\ Controls" "vfe Vote FE"   ///
						"sfe Number of Sons FE" "kfe Number of Children FE"  /// 
						"clust1 Legislators" ///
						"clust2 Votes" ///
						"ymean Mean Dep. Var." ///
						"N Observations")  ///
			mgroups("Hawkish Vote (Narrative Method)" "Hawkish Vote (Extrapolation method)" , pattern(1 0 0 0 1 0 0 0 )  ///
			prefix(\multicolumn{@span}{c}{) suffix(})  ///
			span erepeat(\cmidrule(lr){@span}))  ///
			sfmt(%~12s %~12s  %~12s %~12s %~12s  %~12s %10.2f %10.0f  )  /// 
			label obslast nomtitles booktabs
	


*		B. Excluding Main Votes

		local X terms age age2 i.party_vv house0 
		
			reghdfe w_h1d0_mode   	sons_any_draft_arw_b4 ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 1 & vote_pro == ., a(kids_num_now  sons_num_now ) cluster(id_vv voteFE)
							
				estadd local lfe "No"
				estadd local vfe "No"
				estadd local sfe "Yes"
				estadd local kfe "Yes"
				estadd local con "No"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto reg211
	 
	 
			reghdfe w_h1d0_mode   	sons_any_draft_arw_b4 `X' ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 1 & vote_pro == ., a(kids_num_now  sons_num_now ) cluster(id_vv voteFE)
			
				estadd local lfe "No"
				estadd local vfe "No"
				estadd local sfe "Yes"
				estadd local kfe "Yes"
				estadd local con "Yes"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto reg212
				
				
			reghdfe w_h1d0_mode  	sons_any_draft_arw_b4 ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 1 & vote_pro == ., a(kids_num_now  sons_num_now voteFE) cluster(id_vv voteFE)
							
				estadd local lfe "No"
				estadd local vfe "Yes"
				estadd local sfe "Yes"
				estadd local kfe "Yes"
				estadd local con "No"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto reg213
	 
	 
			reghdfe w_h1d0_mode  	sons_any_draft_arw_b4 `X' ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 1 & vote_pro == ., a(kids_num_now  sons_num_now voteFE ) cluster(id_vv voteFE)
			
				estadd local lfe "No"
				estadd local vfe "Yes"
				estadd local sfe "Yes"
				estadd local kfe "Yes"
				estadd local con "Yes"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto reg214
		

			reghdfe w_h1d0_extrap_mode   	sons_any_draft_arw_b4 ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 1 & vote_pro == ., a(kids_num_now  sons_num_now ) cluster(id_vv voteFE)
							
				estadd local lfe "No"
				estadd local vfe "No"
				estadd local sfe "Yes"
				estadd local kfe "Yes"
				estadd local con "No"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto reg221
	 
	 
			reghdfe w_h1d0_extrap_mode   	sons_any_draft_arw_b4 `X' ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 1 & vote_pro == ., a(kids_num_now  sons_num_now ) cluster(id_vv voteFE)
			
				estadd local lfe "No"
				estadd local vfe "No"
				estadd local sfe "Yes"
				estadd local kfe "Yes"
				estadd local con "Yes"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto reg222
				
				
			reghdfe w_h1d0_extrap_mode  	sons_any_draft_arw_b4 ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 1 & vote_pro == ., a(kids_num_now  sons_num_now voteFE) cluster(id_vv voteFE)
							
				estadd local lfe "No"
				estadd local vfe "Yes"
				estadd local sfe "Yes"
				estadd local kfe "Yes"
				estadd local con "No"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto reg223
	 
	 
			reghdfe w_h1d0_extrap_mode  	sons_any_draft_arw_b4 `X' ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 1 & vote_pro == ., a(kids_num_now  sons_num_now voteFE ) cluster(id_vv voteFE)
			
				estadd local lfe "No"
				estadd local vfe "Yes"
				estadd local sfe "Yes"
				estadd local kfe "Yes"
				estadd local con "Yes"
				estadd sca clust1 = e(N_clust1)
				estadd sca clust2 = e(N_clust2)
				estadd ysumm
				
				est sto reg224
		
		
			esttab   reg211 reg212 reg213 reg214 reg221 reg222 reg223 reg224 using "$output/JPE_REP_TA9b.tex", replace se  drop(terms age age2 terms age age2 *part* *cons* *house*  ) ///
			order("\underline{\emph{All Draft Votes Excluding Main Votes}} \vspace{-0.4cm}" *) ///
			star(* 0.100 ** 0.050 *** 0.010) b(%10.4f %10.4f)  nonotes   ///
			scalars(	"con \\ Controls" ///
					"vfe Vote FE"   ///
					"sfe Number of Sons FE" ///
					"kfe Number of Children FE"  /// 
					"clust1 Legislators" ///
					"clust2 Votes" ///
					"ymean Mean Dep. Var." ///
					"N Observations")  ///
			mgroups("Hawkish Vote (Narrative Method)" "Hawkish Vote (Extrapolation Method)" , pattern(1 0 0 0 1 0 0 0 )  ///
			prefix(\multicolumn{@span}{c}{) suffix(})  ///
			span erepeat(\cmidrule(lr){@span}))  ///
			sfmt( %~12s %~12s  %~12s %~12s %~12s  %~12s  %10.2f %10.0f  )  /// 
			label obslast nomtitles booktabs
		
	
*		C. Non-Draft Votes

		local X terms age age2 i.party_vv house0 

			reghdfe w_h1d0_mode   	sons_any_draft_arw_b4 ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 0, a(kids_num_now  sons_num_now ) cluster(id_vv voteFE)
						
			estadd local lfe "No"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg311
 
 
			reghdfe w_h1d0_mode   	sons_any_draft_arw_b4 `X' ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 0, a(kids_num_now  sons_num_now ) cluster(id_vv voteFE)
 		
			estadd local lfe "No"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg312
			
			
			reghdfe w_h1d0_mode  	sons_any_draft_arw_b4 ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 0, a(kids_num_now  sons_num_now voteFE) cluster(id_vv voteFE)
						
			estadd local lfe "No"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg313
 
 
			reghdfe w_h1d0_mode  	sons_any_draft_arw_b4 `X' ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 0, a(kids_num_now  sons_num_now voteFE ) cluster(id_vv voteFE)
 		
			estadd local lfe "No"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg314
	

			reghdfe w_h1d0_extrap_mode   	sons_any_draft_arw_b4 ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 0, a(kids_num_now  sons_num_now ) cluster(id_vv voteFE)
						
			estadd local lfe "No"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg321
 
 
			reghdfe w_h1d0_extrap_mode   	sons_any_draft_arw_b4 `X' ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 0, a(kids_num_now  sons_num_now ) cluster(id_vv voteFE)
 		
			estadd local lfe "No"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg322
			
			
			reghdfe w_h1d0_extrap_mode  	sons_any_draft_arw_b4 ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 0, a(kids_num_now  sons_num_now voteFE) cluster(id_vv voteFE)
						
			estadd local lfe "No"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg323
 
 
			reghdfe w_h1d0_extrap_mode  	sons_any_draft_arw_b4 `X' ///
				if kids_any_draft_arw_b4 == 1 &  all_draft == 0, a(kids_num_now  sons_num_now voteFE ) cluster(id_vv voteFE)
 		
			estadd local lfe "No"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg324
	
	
			esttab  reg311 reg312 reg313 reg314 reg321 reg322 reg323 reg324 using "$output/JPE_REP_TA9c.tex", replace se  drop(terms age age2 terms age age2 *part* *cons* *house*  ) ///
			order("\underline{\emph{Non-Draft Votes}} \vspace{-0.4cm}" *) ///
			star(* 0.100 ** 0.050 *** 0.010) b(%10.4f %10.4f)  nonotes   ///
			scalars(	"con \\ Controls" ///
						"vfe Vote FE"   ///
						"sfe Number of Sons FE" ///
						"kfe Number of Children FE"  /// 
						"clust1 Legislators" ///
						"clust2 Votes" ///
						"ymean Mean Dep. Var." ///
						"N Observations")  ///
			mgroups("Hawkish Vote (Narrative Method)" "Hawkish Vote (Extrapolation Method)" , pattern(1 0 0 0 1 0 0 0 )  ///
			prefix(\multicolumn{@span}{c}{) suffix(})  ///
			span erepeat(\cmidrule(lr){@span}))  ///
			sfmt( %~12s %~12s  %~12s %~12s %~12s  %~12s %10.2f %10.0f  )  /// 
			label obslast nomtitles booktabs
	
		
*--------------
*		FIG A6
*--------------

*	A. Narrative
		coefplot 	reg111 reg112 reg113 reg114, keep(sons_any_draft_arw_b4)  bylabel("All Draft Votes")  || ///
					reg211 reg212 reg213 reg214, keep(sons_any_draft_arw_b4)  bylabel("Draft Votes Excluding Main Votes") || ///
					reg311 reg312 reg313 reg314, keep(sons_any_draft_arw_b4)  bylabel("Non-Draft Votes") ||,  ///
		byopts(compact cols(1)) xline(0) subtitle(, size(medium) margin(small) justification(left) color(black) bcolor(white) bmargin(top_bottom)) ///
		legend(rows(1)) xtitle("Effect of Draft Age Son on Hawkish Vote (Narrative Method)")  levels(90)  xline(0) xlabel(-0.2(0.1)0.2) msymbol(o)

		/* Save Narrative
		graph save  "$output/JPE_REP_FA6a.gph", replace
		graph export "$output/JPE_REP_FA6a.eps", as(eps) preview(off) replace
		graph export "$output/JPE_REP_FA6a.png", as(png) replace
		graph export "$output/JPE_REP_FA6a.pdf", as(pdf) replace
		*/


*	B. Extrapolation
		coefplot 	reg121 reg122 reg123 reg124, keep(sons_any_draft_arw_b4)  bylabel("All Draft Votes")  || ///
					reg221 reg222 reg223 reg224, keep(sons_any_draft_arw_b4)  bylabel("Draft Votes Excluding Main Votes") || ///
					reg321 reg322 reg323 reg324, keep(sons_any_draft_arw_b4)  bylabel("Non-Draft Votes") ||,  ///
		byopts(compact cols(1)) xline(0) subtitle(, size(medium) margin(small) justification(left) color(black) bcolor(white) bmargin(top_bottom)) ///
		legend(rows(1)) xtitle("Effect of Draft Age Son on Hawkish Vote (Extrapolation Method)")  levels(90)  xline(0) xlabel(-0.2(0.1)0.2) msymbol(o)

		/* Save Extrapolation
		graph save  "$output/JPE_REP_FA6b.gph", replace
		graph export "$output/JPE_REP_FA6b.eps", as(eps) preview(off) replace
		graph export "$output/JPE_REP_FA6b.png", as(png) replace
		graph export "$output/JPE_REP_FA6b.pdf", as(pdf) replace
		*/

	
*--------------
*		TABLE A12
*--------------	

*		A. All Draft Votes

		est clear
		
		local X terms age age2 i.party_vv house0 
		
			reghdfe w_h1d0_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 1, a(kids_num_now  sons_num_now id_vv) cluster(id_vv voteFE)
						
			estadd local lfe "Yes"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg111
 
 
			reghdfe w_h1d0_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 `X' ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 1, a(kids_num_now  sons_num_now id_vv) cluster(id_vv voteFE)
 		
			estadd local lfe "Yes"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg112

			
			reghdfe w_h1d0_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 1, a(kids_num_now  sons_num_now id_vv voteFE) cluster(id_vv voteFE)
						
			estadd local lfe "Yes"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg113
 
 
			reghdfe w_h1d0_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 `X' ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 1, a(kids_num_now  sons_num_now id_vv voteFE) cluster(id_vv voteFE)
 		
			estadd local lfe "Yes"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg114
			

			reghdfe w_h1d0_extrap_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 1, a(kids_num_now  sons_num_now id_vv) cluster(id_vv voteFE)
						
			estadd local lfe "Yes"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg121
 
 
			reghdfe w_h1d0_extrap_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 `X' ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 1, a(kids_num_now  sons_num_now id_vv) cluster(id_vv voteFE)
 		
			estadd local lfe "Yes"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg122

			
			reghdfe w_h1d0_extrap_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 1, a(kids_num_now  sons_num_now id_vv voteFE) cluster(id_vv voteFE)
						
			estadd local lfe "Yes"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg123
 
 
			reghdfe w_h1d0_extrap_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 `X' ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 1, a(kids_num_now  sons_num_now id_vv voteFE) cluster(id_vv voteFE)
 		
			estadd local lfe "Yes"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg124
			

			esttab reg111 reg112 reg113 reg114 reg121 reg122 reg123 reg124 using "$output/JPE_REP_TA12a.tex", replace se  drop(terms age age2 terms age age2 *part* *cons* *house*  ) ///
			order("\underline{\emph{All Draft Votes}} \vspace{-0.4cm}" *) ///
			star(* 0.100 ** 0.050 *** 0.010) b(%10.4f %10.4f)  nonotes   ///
			scalars(	"con \\ Controls" "vfe Vote FE" ///
						"lfe Legislator FE"  ///
						"sfe Number of Sons FE" ///
						"kfe Number of Children FE"  /// 
						"clust1 Legislators" ///
						"clust2 Votes" ///
						"ymean Mean Dep. Var." ///
						"N Observations")  ///
			mgroups("Hawkish Vote (Narrative Method)" "Hawkish Vote (Extrapolation method)" , pattern(1 0 0 0 1 0 0 0 )  ///
			prefix(\multicolumn{@span}{c}{) suffix(})  ///
			span erepeat(\cmidrule(lr){@span}))  ///
			sfmt(  %~12s %~12s  %~12s %~12s %~12s %~12s %~12s  %10.2f %10.0f  )  /// 
			label obslast nomtitles booktabs
			

*		B. Excluding Main Votes

		local X terms age age2 i.party_vv house0 

			reghdfe w_h1d0_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 1 & vote_pro == ., a(kids_num_now  sons_num_now id_vv) cluster(id_vv voteFE)
						
			estadd local lfe "Yes"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg211
 
 
			reghdfe w_h1d0_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 `X' ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 1 & vote_pro == ., a(kids_num_now  sons_num_now id_vv) cluster(id_vv voteFE)
 		
			estadd local lfe "Yes"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg212

			
			reghdfe w_h1d0_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 1 & vote_pro == ., a(kids_num_now  sons_num_now id_vv voteFE) cluster(id_vv voteFE)
						
			estadd local lfe "Yes"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg213
 
 
			reghdfe w_h1d0_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 `X' ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 1 & vote_pro == ., a(kids_num_now  sons_num_now id_vv voteFE) cluster(id_vv voteFE)
 		
			estadd local lfe "Yes"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg214
			

			reghdfe w_h1d0_extrap_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 1 & vote_pro == ., a(kids_num_now  sons_num_now id_vv) cluster(id_vv voteFE)
						
			estadd local lfe "Yes"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg221
 
 
			reghdfe w_h1d0_extrap_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 `X' ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 1 & vote_pro == ., a(kids_num_now  sons_num_now id_vv) cluster(id_vv voteFE)
 		
			estadd local lfe "Yes"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg222

			
			reghdfe w_h1d0_extrap_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 1 & vote_pro == ., a(kids_num_now  sons_num_now id_vv voteFE) cluster(id_vv voteFE)
						
			estadd local lfe "Yes"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg223
 
 
			reghdfe w_h1d0_extrap_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 `X' ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 1 & vote_pro == ., a(kids_num_now  sons_num_now id_vv voteFE) cluster(id_vv voteFE)
 		
			estadd local lfe "Yes"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg224
			

			esttab reg211 reg212 reg213 reg214 reg221 reg222 reg223 reg224 using "$output/JPE_REP_TA12b.tex", replace se  drop(terms age age2 terms age age2 *part* *cons* *house*  ) ///
			order("\underline{\emph{All Draft Votes Excluding Main Votes}} \vspace{-0.4cm}" *) ///
			star(* 0.100 ** 0.050 *** 0.010) b(%10.4f %10.4f)  nonotes   ///
			scalars(	"con \\ Controls" ///
						"vfe Vote FE" ///
						"lfe Legislator FE"  ///
						"sfe Number of Sons FE" ///
						"kfe Number of Children FE"  /// 
						"clust1 Legislators" ///
						"clust2 Votes" ///
						"ymean Mean Dep. Var." ///
						"N Observations")  ///
			mgroups("Hawkish Vote (Narrative Method)" "Hawkish Vote (Extrapolation method)" , pattern(1 0 0 0 1 0 0 0 )  ///
			prefix(\multicolumn{@span}{c}{) suffix(})  ///
			span erepeat(\cmidrule(lr){@span}))  ///
			sfmt(%~12s %~12s  %~12s %~12s %~12s %~12s %~12s %10.2f %10.0f  )  /// 
			label obslast nomtitles booktabs
			

*		C. Non-Draft Votes

		local X terms age age2 i.party_vv house0 
		
			reghdfe w_h1d0_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 0, a(kids_num_now  sons_num_now id_vv) cluster(id_vv voteFE)
						
			estadd local lfe "Yes"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg311
 
 
			reghdfe w_h1d0_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 `X' ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 0, a(kids_num_now  sons_num_now id_vv) cluster(id_vv voteFE)
 		
			estadd local lfe "Yes"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg312

			
			reghdfe w_h1d0_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 0, a(kids_num_now  sons_num_now id_vv voteFE) cluster(id_vv voteFE)
						
			estadd local lfe "Yes"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg313
 
 
			reghdfe w_h1d0_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 `X' ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 0, a(kids_num_now  sons_num_now id_vv voteFE) cluster(id_vv voteFE)
 		
			estadd local lfe "Yes"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg314
			

			reghdfe w_h1d0_extrap_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 0, a(kids_num_now  sons_num_now id_vv) cluster(id_vv voteFE)
						
			estadd local lfe "Yes"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg321
 
 
			reghdfe w_h1d0_extrap_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 `X' ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 0, a(kids_num_now  sons_num_now id_vv) cluster(id_vv voteFE)
 		
			estadd local lfe "Yes"
			estadd local vfe "No"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg322

			
			reghdfe w_h1d0_extrap_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 0, a(kids_num_now  sons_num_now id_vv voteFE) cluster(id_vv voteFE)
						
			estadd local lfe "Yes"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "No"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg323
 
 
			reghdfe w_h1d0_extrap_mode  	sons_any_draft_arw_b4 kids_any_draft_arw_b4 `X' ///
				if (HDeverST ==1 | HDeverSC ==1) & all_draft == 0, a(kids_num_now  sons_num_now id_vv voteFE) cluster(id_vv voteFE)
 		
			estadd local lfe "Yes"
			estadd local vfe "Yes"
			estadd local sfe "Yes"
			estadd local kfe "Yes"
			estadd local con "Yes"
			estadd sca clust1 = e(N_clust1)
			estadd sca clust2 = e(N_clust2)
			estadd ysumm
			
			est sto reg324
			
			esttab reg311 reg312 reg313 reg314 reg321 reg322 reg323 reg324 using "$output/JPE_REP_TA12c.tex", replace se  drop(terms age age2 terms age age2 *part* *cons* *house*  ) ///
			order("\underline{\emph{Non-Draft Votes}} \vspace{-0.4cm}" *) ///
			star(* 0.100 ** 0.050 *** 0.010) b(%10.4f %10.4f)  nonotes   ///
			scalars(	"con \\ Controls" "vfe Vote FE" ///
						"lfe Legislator FE"  ///
						"sfe Number of Sons FE" ///
						"kfe Number of Children FE"  /// 
						"clust1 Legislators" ///
						"clust2 Votes" ///
						"ymean Mean Dep. Var." ///
						"N Observations")  ///
			mgroups("Hawkish Vote (Narrative Method)" "Hawkish Vote (Extrapolation method)" , pattern(1 0 0 0 1 0 0 0 )  ///
			prefix(\multicolumn{@span}{c}{) suffix(})  ///
			span erepeat(\cmidrule(lr){@span}))  ///
			sfmt(%~12s %~12s  %~12s %~12s %~12s %~12s %~12s %10.2f %10.0f  )  /// 
			label obslast nomtitles booktabs
		

*--------------
*		FIG A8
*--------------	
	
*	A. Narrative
		coefplot 	reg111 reg112 reg113 reg114, keep(sons_any_draft_arw_b4)  bylabel("All Draft Votes")  || ///
					reg211 reg212 reg213 reg214, keep(sons_any_draft_arw_b4)  bylabel("Draft Votes Excluding Main Votes") || ///
					reg311 reg312 reg313 reg314, keep(sons_any_draft_arw_b4)  bylabel("Non-Draft Votes") ||,  ///
					byopts(compact cols(1)) xline(0) subtitle(, size(medium) margin(small) justification(left) color(black) bcolor(white) bmargin(top_bottom)) ///
					legend(rows(1)) xtitle("Effect of Draft Age Son on Hawkish Vote (Narrative Method)")  levels(90)  xline(0) xlabel(-0.2(0.1)0.2) msymbol(o)

		/* Save Narrative
		graph save  "$output/JPE_REP_FA8a.gph", replace
		graph export "$output/JPE_REP_FA8a.eps", as(eps) preview(off) replace
		graph export "$output/JPE_REP_FA8a.png", as(png) replace
		graph export "$output/JPE_REP_FA8a.pdf", as(pdf) replace
		*/


*	B. Extrapolation
		coefplot 	reg121 reg122 reg123 reg124, keep(sons_any_draft_arw_b4)  bylabel("All Draft Votes")  || ///
					reg221 reg222 reg223 reg224, keep(sons_any_draft_arw_b4)  bylabel("Draft Votes Excluding Main Votes") || ///
					reg321 reg322 reg323 reg324, keep(sons_any_draft_arw_b4)  bylabel("Non-Draft Votes") ||,  ///
					byopts(compact cols(1)) xline(0) subtitle(, size(medium) margin(small) justification(left) color(black) bcolor(white) bmargin(top_bottom)) ///
					legend(rows(1)) xtitle("Effect of Draft Age Son on Hawkish Vote (Extrapolation Method)")  levels(90)  xline(0) xlabel(-0.2(0.1)0.2) msymbol(o)

		/* Save Extrapolation
		graph save  "$output/JPE_REP_FA8b.gph", replace
		graph export "$output/JPE_REP_FA8b.eps", as(eps) preview(off) replace
		graph export "$output/JPE_REP_FA8b.png", as(png) replace
		graph export "$output/JPE_REP_FA8b.pdf", as(pdf) replace
		*/
	
	
*--------------
*		TABLE A15
*--------------	

		est clear
		local X terms age age2 i.party_vv house0 

		foreach i in 1 2 5 {
									
			reghdfe w_h1d0_mode  	electionyear `X' ///
				if warFE`i' == 1 , a(voteFE ) cluster(id_vv voteFE)	
					
					estadd local lfe "No"
					estadd local vfe "Yes"
					estadd local sfe "No"
					estadd local kfe "No"
					estadd local con "Yes"
					estadd sca clust1 = e(N_clust1)
					estadd sca clust2 = e(N_clust2)
					estadd ysumm
					
					est sto reg`i'1
			 
					
			reghdfe w_h1d0_mode  	electionyear `X' ///
				if warFE`i' == 1 , a(id_vv  voteFE ) cluster(id_vv voteFE)	
								
					estadd local lfe "Yes"
					estadd local vfe "Yes"
					estadd local sfe "No"
					estadd local kfe "No"
					estadd local con "Yes"
					estadd sca clust1 = e(N_clust1)
					estadd sca clust2 = e(N_clust2)
					estadd ysumm
					
					est sto reg`i'2
								
			}

		esttab   reg11 reg12  reg21 reg22  reg51 reg52   using "$output/JPE_REP_TA15.tex", replace se noconstant keep(electionyear ) order(electionyear  ) ///
		star(* 0.100 ** 0.050 *** 0.010) b(%10.4f %10.4f)  nonotes   ///
		scalars(	"lfe \\ Legislator FE" ///
					"vfe  Vote FE"  ///
					"sfe Number of Sons FE" ///
					"kfe Number of Children FE" ///
					"con Other Controls" /// 
					"clust1 Legislators"  ///
					"clust2 Votes" ///
					"ymean Mean Dep. Var." ///
					"N Observations")  ///
		mgroups("WW1" "WW2" "Cold War" , pattern(1  0 1 0 1 0)  ///
		prefix(\multicolumn{@span}{c}{) suffix(})  ///
		span erepeat(\cmidrule(lr){@span}))  ///
		sfmt(%~12s %~12s  %~12s %~12s %~12s  %~12s %~12s %10.3f %10.0f  )  /// 
		label obslast nomtitles booktabs

*--------------
*		TABLE A16
*--------------	

		est clear
		local X terms age age2 i.party_vv house0 

		foreach i in 1 2 5 {
									
			reghdfe w_h1d0_extrap_mode  	electionyear `X' ///
				if warFE`i' == 1 , a(voteFE ) cluster(id_vv voteFE)	
					
					estadd local lfe "No"
					estadd local vfe "Yes"
					estadd local sfe "No"
					estadd local kfe "No"
					estadd local con "Yes"
					estadd sca clust1 = e(N_clust1)
					estadd sca clust2 = e(N_clust2)
					estadd ysumm
					
					est sto reg`i'1
			 
					
			reghdfe w_h1d0_extrap_mode  	electionyear `X' ///
				if warFE`i' == 1 , a(id_vv  voteFE ) cluster(id_vv voteFE)	
								
					estadd local lfe "Yes"
					estadd local vfe "Yes"
					estadd local sfe "No"
					estadd local kfe "No"
					estadd local con "Yes"
					estadd sca clust1 = e(N_clust1)
					estadd sca clust2 = e(N_clust2)
					estadd ysumm
					
					est sto reg`i'2
								
			}

		esttab   reg11 reg12  reg21 reg22  reg51 reg52   using "$output/JPE_REP_TA16.tex", replace se noconstant keep(electionyear ) ///
		order(electionyear  ) star(* 0.100 ** 0.050 *** 0.010) b(%10.4f %10.4f)  nonotes   ///
		scalars(	"lfe \\ Legislator FE" ///
					"vfe  Vote FE"  ///
					"sfe Number of Sons FE" ///
					"kfe Number of Children FE" ///
					"con Other Controls" /// 
					"clust1 Legislators"   ///
					"clust2 Votes" ///
					"ymean Mean Dep. Var." ///
					"N Observations")  ///
		mgroups("WW1" "WW2" "Cold War" , pattern(1  0 1 0 1 0)  ///
		prefix(\multicolumn{@span}{c}{) suffix(})  ///
		span erepeat(\cmidrule(lr){@span}))  ///
		sfmt(%~12s %~12s  %~12s %~12s %~12s  %~12s %~12s %10.3f %10.0f  )  /// 
		label obslast nomtitles booktabs

		
*--------------
*		END LOG
*--------------	
		log close
